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Observations of the Hubble Deep Field South with the 
Infrared Space Observatory - II. Associations and star 
formation rates 
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ABSTRACT 

We present results from a deep mid-infrared survey of the Hubble Deep 
k>( ^ Field South (HDF-S) region performed at 7 and 15/im with the CAM 

, instrument on board the Infrared Space Observatory {ISO). We found re- 

■ liable optical/near ~IR associations for 32 of the 35 sources detected in this 

field by Oliver et al. (2002, Paper I): eight of them were identified as stars, 
one is definitely an AGN, a second seems likely to be an AGN, too, while 
the remaining 22 appear to be normal spiral or starburst galaxies. Using 
model spectral energy distributions (SEDs) of similar galaxies, we compare 
methods for estimating the star formation rates (SFRs) in these objects, 
finding that an estimator based on integrated (3~1000/im) infrared lumi- 
nosity reproduces the model SFRs best. Applying this estimator to model 
fits to the SEDs of our 22 spiral and starburst galaxies, we find that they 
are forming stars at rates of ~ 1 — 100 Mq yr~^, with a median value of 
~ 40 Mq yr""'^, assuming an Einstein - de Sitter universe with a Hubble 
constant of 50 km Mpc^^, and star formation taking place according 
to a Salpeter (1955) IMF across the mass range 0.1 — IOOM0. We split 
the redshift range 0.0 < z < 0.6 into two equal-volume bins to compute 
raw estimates of the star formation rate density, , contributed by these 
sources, assuming the same cosmology and IMF as above and computing 
errors based on estimated uncertainties in the SFRs of individual galaxies. 
We compare these results with other estimates of made with the same 
assumptions, showing them to be consistent with the results of Flores et 
al. (1999) from their ISO survey of the CFRS 1415+52 field. However, 
the relatively small volume of our survey means that our estimates 
suffer from a large sampling variance, implying that our results, by them- 
selves, do not place tight constraints on the global mean star formation 
© 2002 RAS + J V 

^ rate density. 
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1 INTRODUCTION 

In an accompanying paper (Oliver et al. 2002, Pa- 
per I) we described the survey we made of the Hub- 
ble Deep Field South (HDF-SQ) region using the 
LW2 (centred at 6.7/im) and LW3 (15^im) fihers 
of the CAM (Cesarsky et al. 1996) instrument on 
board the Infrared Space Observatory (ISO) (Kessler 
et al. 1996) . A prime motivation for this project came 
from the results of previous ISO surveys, such as our 
own survey (Serjeant et al. 1997, Goldschmidt et al. 
1997, Oliver et al. 1997, Mann et al. 1997, Rowan- 
Robinson et al. 1997) of the northern Hubble Deep 
Field (HDF-N, Williams et al. 1996) and that of Flo- 
res et al. (1999) in the Canada-France Redshift Sur- 
vey 1415-1-52 field (Lilly et al. 1995), which revealed 
an infrared luminosity density at 0.5 z ^ 1 im- 
plying a higher star formation rate density at those 
redshifts than indicated by previous optical studies 
(e.g. Lilly et al. 1996, Madau et al. 1996, Connolly et 
al. 1997). 

These studies were facilitated by the availabil- 
ity of multi-wavelength datasets in those well-studied 
fields, which made possible the association of ISO 
sources with galaxies for which redshifts had been de- 
termined through optical spectroscopy or for which 
they could be estimated on the basis of multi-band 
optical/near-infrared photometry. The same factors 
favour the study of the HDF-S region, and in this pa- 
per we discuss the association of our ISO sources in 
the HDF-S with objects in optical and near-infrared 
catalogues of that field, and the estimation of their 
star formation rates. This process is significantly eas- 
ier than was the case for our survey of the northern 
HDF-N, since, as emphasised in Paper I, our ISO 
dataset in the HDF-S is greatly superior, as a re- 
sult of an observational strategy and data reduction 
procedure both revised significantly in the light of de- 
veloping knowledge of the characteristics of the CAM 
instrument and of the source population it probed. 

The plan of the rest of this paper is as follows. In 
Section 0, we describe the existing optical and near- 
infrared catalogues in the HDF-S region with which 
we shall seek associations for our ISO sources, review 
the likelihood ratio method used to obtain them and 
present the results of its application. Section ^ com- 
pares a variety of commonly used star formation rate 
(SFR) estimators, through applying them to model 
SEDs, while, in Section W, we apply the best of them 
to the multi-wavelength photometric datasets com- 
piled for the galaxies with which we have associated 
our 750 sources, to yield estimates of their SFRs, and 
Section ^ uses these results to assess the contribution 
of /50-selected sources to the star formation history 
of the HDF-S. Finally, in Section |^ we discuss the re- 
sults of this paper and the conclusions we draw from 
them. Appendices present mathematical details of the 
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Figure 1. This figure shows the location of our ISO rasters 
with respect to those of other datasets taken in the area. 
The shaded regions mark the HST fields, with the STIS 
and NICMOS fields to the east and south of the WFPC2 
field, respectively. The thick-dashed irregular shape and 
the solid circle show, respectively, the maximum extent of 
our ISO coverage (the coverage in the two bands difl'ers 
slightly) and the region from which the source catalogues 
of Paper I were selected. The remaining lines show bound- 
aries of four optical/near-IR surveys discussed in Section 
^, as follows: (i) dotted line — AAT prime focus imaging 
survey of Verma et al. (2002); (ii) dashed line - CTIO 
BTC survey of Gardner e t al. (1999): ( iii) dot-dashed line 
- CTIO BTC survey of ([Walker 1999); (iv) dot-dot-dot- 



dashed line - ESO EIS optical imaging survey of Da Costa 
et al. (1998); and (v) long dashed line - ESO EIS near- 
infrared survey of da Costa et al. (1998). 



processes of correcting an SFR estimate to a canoni- 
cal IMF, and of computing the sampling variance in 
estimates of the SFR density using a lognormal model 
for the cosmological density field. 



2 OPTICAL CATALOGUES AND 
ASSOCIATIONS 

As shown in Figure ^ the region within which the fi- 
nal source catalogue of Paper I was selected is covered 
by a number of optical and near-infrared imaging sur- 
veys and we have sought to associate our ISO sources 
with objects detected in all of them, using the likeli- 
hood ratio method detailed by Sutherland & Saunders 
(1992). 



2.1 Method 

The implementation of the likelihood ratio associa- 
tion technique used here is essentially identical to that 
we used for the /S'O-HDF-N sources in Mann et al. 
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(1997), so we only briefly review it here. The hke- 
lihood ratio, LR, for the association of a particular 
source with a given catalogue object is defined to be 
the ratio, ptrua/Pchance, of the infinitesimal probabil- 
ity of finding the true counterpart to the source at 
the position of the object and with its flux to the 
inflnitesimal probability of an object with that flux 
being found there by chance. Sutherland & Saunders 
(1992) show that it takes the form 

LRif,.,y)^^-in-^, (1) 

where e{x,y) is the probability distribution for posi- 
tional offsets, {x,y), between source and object [nor- 
malized so that J e(x, y) dx Ay — 1 with the integral 
being taken over all space], n(/) is the surface density 
of objects per unit interval in flux, /, and g(/) is the 
probability distribution function for an ensemble of 
sources, measured in the same passband in which the 
object catalogue is deflned. 

The quantity g(/) essentially measures the flux 
distribution of the true counterparts of the sources in 
the object catalogue, and is unknown unless associa- 
tions have previously been found between an ensem- 
ble of similar sources and a catalogue of similar ob- 
jects. Sutherland & Saunders (1992) suggested that 
it may be estimated in binned form by subtracting 
from the flux histogram of the objects lying near the 
source positions that for the full object catalogue. If 
the latter is scaled to cover the same area as used 
to compute the former, then the resulting quantity is 
proportional to q{f), exhibiting an excess of galax- 
ies in the bins corresponding to the fluxes of the true 
counterparts of the sources in the object catalogues. 
Mann et al. (1997) showed that, when applied to rel- 
atively small source samples, such as in our 750 sur- 
veys of the two Hubble Deep Fields, this empirical 
g(/) is very noisy, and they preferred to take q[f) to 
be a constant, independent of flux. Since their Fig- 
ure 1 strongly implied that the ISO-HDF-N sources 
were associated with objects brighter than J8i4 — 23 
in the catalogue of Williams et al. (1996), this as- 
sumption would have led them to underestimate the 
value of ptruc for bright objects, but, in practice, this 
was seen to have no effect on the associations made, 
so we follow the same procedure here. The validity 
of making that same simplifying assumption here is 
illustrated by Fig, ^ where we compare the R band 
magnitude distribution near (deflned to be within 6" 
of) the 30 750 sources from Table 8 of Paper I con- 
tained within the area of the Gardner et al. catalogue, 
and that for the catalogue as a whole. This shows the 
strong excess of bright (R=19-21) objects near ISO 
source positions - a two-sided Kolmogorov-Smirnov 
test gives a probability of P ~ 10"^° that these two 
are drawn from the same parent distribution - which 
gives us confidence both in the reliability of our source 
detections and in our adopting a constant q(f) here. 
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Figure 2. A comparison of (solid line) the R band mag- 
nitude distribution of the objects within circles of radius 
6" around the 30 ISO sources from Table 8 of Paper I 
contained within the area of the Gardner et al. catalogue, 
and (dashed line) that for the catalogue as a whole, show- 
ing the strong excess of bright objects near ISO source 
positions: a two— sided Kolmogorov— Smirnov test gives a 
probabiUty of P ~ 10"^'' that these two are drawn from 
the same parent distribution. 

For each association made between a given 750- 
HDF-N source and a particular object in the 78i4- 
band catalogue of Williams et al. (1996), Mann et 
al. (1997) computed the probability, Pran(78i4), that 
a fictitious source, placed at random in the HDF-N, 
would have a likeliest association with an object in 
the 7gi4-band catalogue producing a likelihood ratio 
at least as high as that for the source and object in 
question. Clearly, Pran combines information as to the 
reliability of the source detection as well as the proba- 
bility that the given object is the correct counterpart 
of the source, but, as discussed below, it provides a 
useful measure of the quality of the associations made 
in this case: a more direct measure would be available 
(Sutherland & Saunders 1992, Rutledge et al. 2000) 
if we could compute our LR values using a good es- 
timate of q{f), but, as mentioned above, the small 
number of sources here leads us to take q{f) to be 
a constant, thereby leaving LR defined only up to a 
multiplicative factor, in which case we cannot employ 
the formalism of Sutherland & Saunders (1992) or 
Rutledge et al. (2000) for quantifying the reliability 
of associations. 

We sought associations for our 75*0 sources with 
objects in the following set of catalogues: (i) the R 
band AAT prime focus catalogue of Verma et al. 
(2002); (u) the 78i4-band STScI HST cataloguej]; (in) 
the I band EIS catalogue of da Costa et al. (1998); 
(iv) the K band EIS catalogue of da Costa et al. 
(1998); and (v) the R band CTIO BTC catalogue 

t available from the STScI HDFS WWW site at 
www. stsci . edu/ftp/science/hdf south/catalogs .html 
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Figure 3. The Pran histogram for the associations made 
between the 35 sources in Table 8 of Paper I and the R 
band catalogue of Gardner et al. (1999), supplemented by 
the ESO EIS K or I band catalogue for those sources falling 
into the masked region of the Gardner et al. image. We 
define good associations to be those with Pran < 0.15, to 
the left of the dashed line. 



(version 2) of Gardner et al. (1999), available from 
http : //hires .gsfc .nasa.gov/~research/hdf s-btc/. 

Inspection of the results of this process revealed that 
almost all of the ISO sources were associated with 
the same star or galaxy in each case, indicating the 
robustness of this method of association. Moreover, 
the associations made using the different catalogues 
in their different bands yielded very similar Pran val- 
ues for the associations between given pairs of source 
and object, giving us confidence that it is a sound 
statistic to use. 



2.2 Results 

As mentioned above, very similar Pran results are ob- 
tained for our sources when the likelihood ratio pro- 
cedure is run against all five optical and near-IR cat- 
alogues listed in the previous subsection. We take our 
principal association to be with the R band CTIO 
BTC catalogue of Gardner et al. (1999), simply be- 
cause it includes the largest number of the objects 
associated with our sources. In Figure ^ we plot the 
histogram of Pran values resulting from the associa- 
tion of the ISO sources in Table 8 of Paper I with ob- 
jects in the CTIO catalogue. A few of our sources are 
clearly associated with objects in regions of the image 
masked out by Gardner et al. (1999) when they con- 
structed their SExtractor (Bertin & Arnouts 1996) 
catalogue from it: for these cases we use the Pran 
value for associated EIS K band or I band object in- 
stead. This histogram is similar to what we should 
expect: most of the ISO sources have been associated 
with objects with low Pran values, indicating a high 
probability that this is a correct identification, while 



either spurious or have optical counterparts too faint 
for detection in the CTIO image, which has a limiting 
magnitude of P ~ 23. The Pran value to be taken as 
the threshold below which associations are defined to 
be good is somewhat arbitrary when the number of 
sources is low, as here: with a much larger number 
it is possible to see the point where the number of 
sources begins significantly to exceed the fiat tail due 
to sources which are spurious or have no counterpart 
in the object catalogue, but that is clearly not possi- 
ble in Figure ^ The form of that histogram suggests 
that Pran = 0.15 might be an appropriate threshold. 
Note that Mann et al. (1997) took their threshold for 
reasonable identifications of ISO-HDF-N sources to 
be 0.35, so the much tighter bunching to lower Pran 
values here indicates once more that the ISO data 
presented here are of a far higher quality than those 
used in the initial ISO-HDF-N analysis. In partic- 
ular we have benefitted from the knowledge of the 
ISO-CKiA image distortion which was uncharacter- 
ized at the time that Mann et al. (1997) performed 
their identification of the optical counterparts of the 
ISO-HDF-N sources. 

In Figure ^ we show postage stamp images for 
the 32 sources in Paper I with reliable associations, 
overlaying ISO 6.7 and 15 /xm contours onto opti- 
cal images, with the associated object lying at (0,0). 
Table hi summarizes the properties of these objects, 
while the next subsection presents more details of the 
associations. As mentioned previously, principal asso- 
ciations were made with the GSFC R band catalogue 
of Gardner et al. (1999), so the optical magnitudes 
tabulated are from their catalogue, with one excep- 
tion (ISOHDFS J223237-603235), which is omitted 
from their catalogue, causing us to use the EIS op- 
tical catalogue of da Costa et al. (1998). The photo- 
metric data in the GSFC catalogue are calibrated us- 
ing standards measured in Johnson UBV and Cousins 
RI, so we quote all optical magnitudes in that sys- 
tem. For the one case of ISOHDFS J223237-603235, 
this involves converting the UBVRI magnitudes back 
to the Johnson- Cousins system from the AB system 
in which da Costa et al. (1998) present their pho- 
tometry, using the conversions given in their paper, 
namely U = Uab - 0.82, B = Bab + 0.06, V = Vab, 
R = Pab - 0.17 and / = Jab - 0.42. For their near- 
infrared data, da Costa et al. (1998) performed photo- 
metric calibration in JHK using Persson et al. (1998) 
standards and then converted them to AB magni- 
tudes: again we have reversed this conversion, using 
the relations quoted by da Costa et al. (1998), namely 
J = Jab -0.89, H = Ji-AB-1.38 and K = A'ab-1.86. 
In Table |l] we present spectroscopic redshifts from 
Glazebrook et al. (2002) where available, and pho- 
tometric redshifts where not. Photometric redshifts 
were estimated by three of us, using independent 
methods: M. Rowan-Robinson (MRR) used (Rowan- 
Robinson 2001) SEDs based on those of Yoshii & 
Takahara (1988), A. Franceschini (AF) used the PE- 
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Figure 4. Postage stamp images for the 32 sources from Paper I witli reliable associations. ISO contours from the 
6.7/xm (dashed) and 15fim (solid) signal/noise map are plotted over an optical/near-infrared image, centred on the object 
associated with the ISO source: the title of each image gives the name of the ISO source and the optical/NIR image used 
for the background. High resolution version of this figure available from astro.ic.ac.uk/hdfs. 



Figure 4 - continued Postage stamp images for the 32 sources from Paper I with reliable associations. ISO contours 
(1, 2, 3 . . . 9, 10, 20, . . . , 90, 100, 200) from the 6.7fim (dashed) and 15/im (solid) signal/noise map are plotted over an 
optical/near- infrared image, centred on the object associated with the ISO source: the title of each image gives the name of 
the ISO source and the optical/NIR image used for the background. High resolution version of this figure available 
from astro.ic.ac.uk/hdfs. 



GASE models of Fioc & Rocca-Volmerange (1997) 
and R.G. Mann (RGM) used the GRASIL (Silva et al. 
1998) model fits to SEDs of a range of known galax- 
ies (Arp220, MlOO, M51, M82, NGC6090, NGC6946) 
presented by Silva et al. (1998). Results for individ- 
ual galaxies are discussed below, but, in most cases, 
these independent methods yielded photometric red- 
shift estimates that agreed with each other, and with 
those computed by Stephen Gwyn (Gwyn 1999) and 

the SUNY group|^, to 5z ~ 0.1, which is therefore 
the accuracy we claim for the adopted photometric 
redshifts listed in parentheses in Table hi. 



2.3 Notes on individual sources 

(i) ISOHDFS J223237-603235: The association 
here is with an 1=16 star, 2 arcsec from the ISO 
source position and yielding Pran = 0.001 from the 
EIS optical catalogue (the star is masked out of the 
CTIO catalogue). 

(ii) ISOHDFS J223237-603256: As shown by 
the contours in Figure ^ this 6.7 /im source is very 
safely (Pran ~ 0.001) associated with a bright (17th 
magnitude in R) object clearly identified as a star by 
SExtractor in the AAT R and EIS I images. 

(iii) ISOHDFS J223240-603141: This 15 ^im 
source is associated fairly securely (Pran = 0.095) with 
an 1=20 galaxy, located 3 arcsec from the ISO posi- 
tion. MRR estimates Zphot = 2.02 for this galaxy, on 
the basis of its UBVRI magnitudes, while both AF 
and RGM estimate 0.45, so it is this latter figure we 
adopt: as shown in Figure ^, with this assumed red- 
shift, the optical and ISO data for this galaxy are well 
fitted by the GRASIL model for NGC6090. 

(iv) ISOHDFS J223243-603242: This is an 
1=20 galaxy, detected significantly (> 4a) in the co- 
added maps at both 6.7 and 15 ^im. This association 
yields Pran = 0.047, and MRR and AF both estimate 
a photometric redshift of 2:phot = 0.5 from the op- 
tical/NIR data for this galaxy, while RGM obtains a 
slightly lower value of 0.45: we adopt 0.5, and show in 
Figure H that with this redshift, the U band to 6.7/im 



i See www.ess.sunysb.edu/astro/hdfs/home.html 
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data are in good agreement with the GRASIL model 
for the starburst galaxies (Arp220, M82, NGC6090). 

(v) ISOHDFS J223243-603351: Another very 
significant detection at both 6.7 and 15 fim, this 
source is associated (Pran = 0.028) with a 19th mag- 
nitude galaxy, with detections at both 4.9 and 8.6 
GHz (with fiuxes of 0.163 and 0.111 mjy respectively: 
A. Hopkins, priv. comm.). Its spectrum exhibits one 
broad line, and yields a redshift of z = 0.0918, so this 
appears to be a low-redshift AGN: RGM, AF, MRR 
estimated Zphot ~ 0.00,0.15,0.12, respectively. 

(vi) ISOHDFS J223243-603441: This source, 
detected very significantly in both ISO bands, is asso- 
ciated (Pran = 0.021) with the brighter (19th magni- 
tude in I) of a close pair of galaxies, for which RGM, 
AF, MRR estimate Zphot = 0.40,0.50,0.59, respec- 
tively, on the basis of UBVRIJHK photometry. We 
adopt z = 0.5 and show in Figure H that this gives 
this galaxy an SED more like the cirrus galaxy M51 
than the starbursts. 

(vn) ISOHDFS J223244-603110: This source, 
solidly detected at 15 /im, is located 3 arcsec from 
an 1=20 galaxy, with which we associate it securely 
(Pran = 0.066). The UBVRIJHK photometry for this 
galaxy yields a fairly wide spread of photometric red- 
shift estimates, with RGM, AF and MRR obtaining 
0.10, 0.25 and 0.32, respectively. We adopt the mid- 
dle value, using which the galaxy's SED is seen to 
be similar to that of the GRASIL fit to the starburst 
NGC6090, and note that this is one of our more un- 
certain redshift estimates. 

(viii) ISOHDFS J223244-603455: This 15^im 
source is associated with an 1=19 galaxy 2 arcsec from 
the ISO position, with a Pran value of 0.031, and UB- 
VRIJHK magnitudes yielding photometric redshift 
estimates of 0.32, 0.30 and 0.45 from MRR, RGM and 
AF, respectively. A value of «phot = 0.35 is adopted, 
with which the SED of this galaxy matches that of 
the GRASIL model for the normal spiral MlOO, as 
shown in Figure ^. 

(ix) ISOHDFS J223245-603226: This 15 fim 
source is associated with one of a pair of close (in- 
teracting?) galaxies, with 1=20 and Pran = 0.080, 
for which RGM, AF and MRR estimate photomet- 
ric redshifts of 0.50, 0.50 and 0.51. Rigopoulou et al. 
(2000) measured a spectroscopic redshift of 2 = 0.59, 
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Figure 4 — continued Postage stamp images for the 32 sources from Paper I with reliable associations. ISO contours 
(1, 2, 3 . . . 9, 10, 20, . . . , 90, 100, 200) from the 6.7/im (dashed) and 15/im (solid) signal/noise map are plotted over an 
optical/near-infrared image, centred on the object associated with the ISO source: the title of each imago gives the name of 
the ISO source and the optical/NIR image used for the background. High resolution version of this figure available 
from astro.ic.ac.uk/hdfs. 



Figure 4 - continued Postage stamp images for the 32 sources from Paper I with reliable associations. ISO contours 
(1, 2, 3 . . . 9, 10, 20, . . . , 90, 100, 200) from the 6.7/im (dashed) and 15/im (solid) signal/noise map are plotted over an 
optical/near-infrared image, centred on the object associated with the ISO source: the title of each image gives the name of 
the ISO source and the optical/NIR image used for the background. High resolution version of this figure available 
from astro.ic.ac.uk/hdfs. 



and, using that value, this galaxy's SED is a reason- 
able fit to the GRASIL model SED for the starburst 
NGC6090. 

(x) ISOHDFS J223245-603418: This is an 1=20 
galaxy, detected at both 6.7 and 15 /im, yielding 
Pjan = 0.076, with a spectrum showing H/3 and 
OllI (4959+5007) lines, from which a redshift of 2 = 
0.4606 was determined: RGM, AF, MRR estimated 
Zphot = 0.55, 0.6, 0.95, respectively, for this source. 
This is also a radio source, detected at 1.4, 2.5 and 
4.9 GHz, with fluxes of 0.200, 0.149 and 0.127 mjy, 
respectively (A. Hopkins, priv. comm.). As shown by 
Figure ^, the SED of this galaxy is in good agreement 
with the GRASIL model for Arp220 over six decades 
in wavelength. 

(xi) ISOHDFS J223247-603335: Detected at 
both 6.7 and 15 /^m, we associate (Pran = 0.029) this 
source with an 1=19 galaxy, with a spectrum display- 
ing a number of emission lines and yielding a redshift 
of 2 = 0.5803: RGM AF, MRR, Gwyn and SUNY 
group estimated photometric redshifts of 0.50, 0.60, 
0.52, 0.56, 0.66, respectively, for this galaxy, whose 
SED is similar to that of the GRASIL model for the 
spiral galaxy M51. 

(xii) ISOHDFS J223250-603359: This 6.7 /im 
is associated very securely (Pran = 0.003) with an 
1=17 M2V star displaying clear diffraction spikes in 
its WFPC2 image. 

(xin) ISOHDFS J223251-603335: This is iden- 
tified with an 1=22 galaxy 3 arcsec away, yielding 
Pran = 0.102 for associatiou with the EIS K band cat- 
alogue. RGM, AF, MRR, Gwyn and the SUNY group 
estimate Zph= 0.50, 0.70, 0.95, 0.56, 0.57, respectively, 
for this galaxy, whose SED (with an adopted redshift 
of 0.7) shows the rise through the infrared character- 
istic of a starburst galaxy. 

(xiv) ISOHDFS J223252-603327: This is a ra- 
dio source (0.109 mJy at 1.4 GHz: A. Hopkins, priv. 
comm.), located at the centre of quite an extended 
region of emission in both ISO bands. The source is 
associated (Pran = 0.073) with the EIS K band cat- 
alogue. SExtractor stellarity indices exist for this ob- 
ject in U, B, V, R, I, J, H and K bands from the EIS 
catalogue, and vary from 0.48 in U to 0.97, indicat- 
ing that the image is getting more point-like at longer 
wavelengths. This, together with the radio detection. 



suggests that this may be an obscured AGN, although 
Rigpoulou et al. (2000) do not report any strong AGN 
features in their spectrum of this object. That yielded 
z = 1.27, in excellent agreement with the photomet- 
ric redshifts of Gwyn and the SUNY group (1.28 and 
1.27, respectively), while RGM, AF and MRR ob- 
tained more widely-varying values of 0.50, 0.90 and 
1.46 for this galaxy. 

(xv) ISOHDFS J223254-603115: This source, 
detected at both 6.7 and 15 /im, is associated with the 
brighter of a pair of very close (interacting?) galaxies 
separated by 2 arcsec. In the ESO EIS near-infrared 
catalogue of da Costa et al. (1998), the two galax- 
ies have (J,H,K) magnitudes of (21.12, 20.64, 20.43) 
and (21.63, 21.18, 21.03), while the brighter galaxy 
has 1=19.71 in the GSFC CTIO catalogue, yielding 
Pran = 0.044, although it appears that the optical 
magnitudes include a significant contribution from 
the second galaxy, and, as a result the optical and 
near-infrared photometry in Figure ^ do not appear 
to be consistent. Fortunately, a spectroscopic redshift 
of z — 0.5111 was measured for this galaxy by Glaze- 
brook et al. (2002), so it is not necessary to try to esti- 
mate a photometric redshift from this inconsistent set 
of magnitudes. Assuming that the near-infrared pho- 
tometry is correct, the SED obtained using z = 0.5111 
is very similar to that for the GRASIL model for the 
starburst galaxy NGC6090. 

(xvi) ISOHDFS J223254-603129: This is an in- 
triguing case. The peaks of 6.7 and 15/im emission are 
very close to each other, but 4-5 arcsec away (in direc- 
tions 120° apart) from both an 1=21 galaxy (yielding 
Pran = 0.162) and a 0.143 mJy 1.4 GHz radio source, 
which appear not to be associated. It is possible that 
the ISO source is not associated with that galaxy ei- 
ther (but, instead, with an object too faint for these 
optical/near-infrared survey data that may or may 
not be the source of the radio emission), but we as- 
sume here that it is. RGM, MRR and AF estimated 
redshifts of 0.15, 0.48 and 0.50 for this galaxy, and, 
adopting Zphot = 0.2 we see, from Figure |^, that we 
obtain an SED similar to that of the GRASIL model 
for Arp220. 

(xvii) ISOHDFS J223254-603143: This 6.7Aim 
source is associated securely (Pran = 0.009) with an 
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Figure 4 - continued Postage stamp images for the 32 sources from Paper I with reliable associations. ISO contours 
(1, 2, 3 . . . 9, 10, 20, . . . , 90, 100, 200) from the d.Tfim (dashed) and 15fim (solid) signal/noise map are plotted over an 
optical/near-infrared image, centred on the object associated with the ISO source: the title of each image gives the name of 
the ISO source and the optical/NIR image used for the background. High resolution version of this figure available 
from astro.ic.ac.uk/hdfs. 



Figure 4 - continued Postage stamp images for the 32 sources from Paper I with reliable associations. ISO contours 
(1, 2, 3 . . . 9, 10, 20, . . . , 90, 100, 200) from the 6.7^m (dotted) and Ibfim (solid) signal/noise map are plotted over an 
optical/near-infrared image, centred on the object associated with the ISO source: the title of each image gives the name of 
the ISO source and the optical/NIR image used for the background. High resolution version of this figure available 
from astro.ic.ac.uk/hdfs. 



1=18 star, which shows clear diffraction spikes in HST 
imaging data. 

(xviii) ISOHDFS J223256-603059: This 6.7 
source has no reasonable association: its best associ- 
ation in the GSFC CTIO catalogue is with an 1=22.5 
galaxy 5 arcsec away, yielding Pran = 0.398. The 
source is located very close to the very bright source 
ISOHDFS J223259-603118, and may be an artifact 
produced by inaccuracies in the application of our 
background subtraction procedure so close to this, the 
second brightest 6.7 fim source in our catalogue. 

(xix) ISOHDFS J223256-603513: There is no 
reasonable association for this 15 /xm source in the 
GSFC CTIO catalogue, which is the only one covering 
this area: the best association is with an 1=22 galaxy 
7 arcsec away, which yields Pran = 0.697. 

(xx) ISOHDFS J223257-603305: This 15 ^im 
source is associated (Pran = 0.047) with an 1=20 
galaxy for which Glazebrook et al. (2002) determined 
a redshift of z = 0.5823 from a spectrum exhibiting a 
number of narrow emission lines, which, together with 
the relative levels of the 15/im detection and 6.7/im 
upper limit, suggest that this is a starburst galaxy, 
although its SED is not a particularly good match to 
any of the GRASIL models. RGM, AF, MRR had es- 
timated Zphot = 0.40, 0.60, 0.62, respectively, for this 
source. 

(xxi) ISOHDFS J223259-603118: This strong 
6.7/im source is a bright star, detected also at 15 /im 
but masked out of most of the SExtractor catalogues 
created from the optical/near-IR surveys under dis- 
cussion here. It appears in the EIS survey as a K=14 
G2III star, yielding Pran < 0.001. 

(xxii) ISOHDFS J223302-603137: Like ISO- 
HDFS J223256-603059, this 6.7 /im source has no 
reasonable identification: the best candidate in the 
GSFC CTIO catalogue is a faint (1=25) galaxy lying 
almost 5 arcsec away, and yielding Pran = 0.539. This 
source lies close to the bright star ISOHDFS J223259- 
603118, and, like the other unidentified 6.7/im source, 
ISOHDFS J223256-603059, may be an artifact resuh- 
ing from background subtraction errors. 

(xxiu) ISOHDFS J223302-603213: The object 
associated with this source is bright (1=17) and less 
than 1 arcsec from the ISO position, yielding a Pran 
value of 0.006. It is classed as stellar by SExtractor, 



and classified as an M3V star on the basis of the spec- 
trum taken by Glazebrook et al. (2002). 

(xxiv) ISOHDFS J223303-603230: The associ- 
ation for this source, detected significantly at both 6.7 
and 15 /im, is a bright object, spectroscopically con- 
firmed to be an MIV star using the Glazebrook et al. 
(2002) spectrum. It is masked out of the GSFC CTIO 
catalogue, but has K=15 in the EIS near-infrared cat- 
alogue, yielding Pran = 0.001. 

(xxv) ISOHDFS J223302-603323: This source, 
detected in both bands, is associated (Pran = 0.049) 
with an 1=20 galaxy for which RGM, AF, MRR, 
Gwyn and the SUNY group estimate photometric 
redshifts of 0.40, 0.60, 0.66, 0.474 and 0.400, respec- 
tively. With an adopted Zphot ~ 0.60, we obtain an 
SED similar to that of the GRASIL models for nor- 
mal spirals like MlOO, NGC6946 or M51. 

(xxvi) ISOHDFS J223303-603336: The 1=20 
galaxy associated with this source is 6 arcsec from 
the ISO position, which is the cause of its relatively 
poor Pran value of 0.131, but the ISO source posi- 
tion could be shifted due to the close proximity of the 
bright source ISOHDFS J223306-603349 (the bright- 
est in our 15/im catalogue): once again, our back- 
ground subtraction method could be leaving artifacts 
close to this bright source. RGM, AF, MRR, Gwyn 
and the SUNY group estimate photometric redshifts 
of 0.30, 0.35, 0.35, 0.419 and 0.440, respectively, for 
this galaxy, and with an adopted redshift of 0.35 we 
obtain an SED which fits the GRASIL starburst mod- 
els well from the U band to 6.7/im: we assume that 
the absence of a detection at 15/im is due to prob- 
lems with the subtraction of the side-lobes of ISO- 
HDFS J223306-603349. 

(xxvii) ISOHDFS J223306-603349: This 1=16 
spiral galaxy, with a spectroscopic redshift deter- 
mined by Glazebrook et al. (2002) to be z = 0.1733, is 
the brightest 15 /xm source in our catalogue, and has 
Pran = 0.002. It is detected in the radio at 0.533 and 
0.300 mjy at 1.4 and 2.5 GHz, respectively (A. Hop- 
kins, priv. comm.), and, from Figure ^ we see that 
this galaxy has the SED of a normal spiral galaxy, 
rather than a starburst. RGM, AF, MRR estimated 
2phot = 0.15,0.25,0.15, respectively, for this galaxy. 

(xxviu) ISOHDFS J223306-603436: This 1=20 
galaxy yields a Pran value of 0.095, and is detected 
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Table 1. Properties of the objects associated with the ISO sources in the Hubble Deep Field South. 
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Notes to Table gj: The positions and Pran values tabulated refer to the associated object in the GSFC catalogue, with the exception 
of the 8 sources which have Pran values in parentheses, which denote that they refer to the association with an object in the EIS 
near-infrared catalogue. For five of these cases (those with positions in parentheses) this is because the source is missing from the 
GSFC catalogue, while the remaining have significantly lower Pran values resulting from their likelihood ratio association with the 
EIS K band catalogue than with the GSFC R band catalogue, as a result of red R — K colours: these three objects are the first, 
second and ninth reddest of the 24 associated objects for which R — K has been measured. 



° All RAs prefixed by 22''. 
' All Decs prefixed by -60°. 

The object associated with this source is not included in the GSFC catalogue, so the tabulated photometric data come from the 
EIS catalogue of da Costa et al. (1998) 

The position of the object associated with this source in the GSFC catalogue suggests that the optical magnitudes quoted here 
may be too bright, through the inclusion of a close, faint companion, which the EIS near-infrared catalogue marks as a separate 
object, but which the GSFC catalogue does not. 



in both 750 bands. MRR, RGM and AF estimate 
photometric redshifts of 0.35, 0.60 and 1.00 for this 
galaxy, and an adopted Zphot ~ 0.60 gives an SED 
very similar to that of the GRASIL models for the 
starburst galaxies M82 and NGC6090. 

(xxix) ISOHDFS J223306-603450: The identi- 
fication of this source is complicated. Its ISO source 
position is less than 1 arcsec from a faint galaxy 
(1=22), which is found with K=21 in the EIS cata- 
logue. The Pran value resulting from the association 
of this object in the GSFC CTIO catalogue is 0.307, 



but it is a much better 0.094 when computed in the 
K band, indicating that this galaxy is very red in 
R-K. RGM, MRR and AF estimate redshifts for it of 
0.65, 1.05 and 0.90, and with an adopted Zphot ~ 0.75 
we obtain an SED similar to the GRASIL model of 
Arp220. 

(xxx) ISOHDFS J223307-603248: This source, 
detected in both bands, is securely (Pran ~ 0.034) as- 
sociated with an 1=20 galaxy for which Glazebrook 
et al. (2002) have determined a spectroscopic red- 
shift of z — 0.513, yielding an SED similar to the 
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GRASIL model for the starburst NGC6090. RGM, 
AF, MRR estimated Zphot = 0.40, 0.50, 0.55, respec- 
tively, for this galaxy. 

(xxxi) ISOHDFS J223308-603314: The likeli- 
hood ratio procedure associates (with P^an ~ 0.034) 
this source with an 1=17 K4V star, although the ISO 
contours in both bands are centred closer to an 1=20 
galaxy 5 arcsec away. AF, MRR and RGM estimate 
redshifts of 0.05, 0.10 and 0.50, and we choose the 
last of these, as that gives the best fit to one of the 
GRASIL SEDs, that for the NGC6090 model, but 
highlight that this is one of our most uncertain red- 
shift estimates. 

(xxxii) ISOHDFS J223312-603350: Detected 
significantly in both bands, this source is associated 
(with Pran = 0.103) with an 1=20 galaxy 2 arcsec 
from the 6.7^m source position and 3 arcsec from 
that at 15pm. Glazebrook et al. (2002) took a spec- 
trum at this sky position, but it yielded no features 
capable of determining the galaxy's redshift. RGM, 
AF and MRR determined photometric redshifts of 
0.45, 0.70 and 0.32, respectively, but none of these 
yields an SED in particularly good agreement with 
the GRASIL starburst models: we adopt Zphot = 0.50, 
as that gives, perhaps, the best agreement (with the 
Arp220 SED), but note that this is highly uncertain. 

(xxxiu) ISOHDFS J223312-603416: This ISO 
source position lies 1 arcsec from a faint (K=21) EIS 
galaxy, which is not present in the GSFC CTIO cat- 
alogue. This near-infrared association yields Pran = 
0.128, and MRR, AF and RGM estimate redshifts 
of 1.24, 1.30 and 2.05, respectively, although this is 
clearly very uncertain, since it is determined from 
JHK photometry alone. The relative fluxes in the two 
ISO and three NIR bands suggests that this is a star- 
burst galaxy, but none of the estimated redshifts listed 
above results in a good fit to any of the GRASIL SED 
models; for definiteness we assume 2;phot = 1-3, but 
stress that this is not at all well-constrained. 

(xxxiv) ISOHDFS J223314-603203: This is an- 
other example of a background subtraction arti- 
fact complicating the identification of a source. This 
source, detected significantly in both bands, lies close 
to a very bright star (the brightest 6.7 fim source 
in our catalogue, and the second brightest at 15/im) 
and has no reliable association. The region is masked 
out of the GSFC CTIO catalogue, while the EIS 
near-infrared catalogue produces an association with 
Pjjjjj = 0.111 with a K=21 galaxy 3 arcsec from the 
ISO source position, but the contours in this region 
have probably been disturbed by the background sub- 
traction procedure. RGM and MRR estimated red- 
shifts of 0.45 and 0.82, respectively, for this galaxy, 
on the basis of JHK photometry only, while we note 
that a significantly lower redshift of ^phot = 0.20 pro- 
duces a reasonable match to the SED of the starburst 
NGC6090, so we adopt that here: it is clear, however, 
that this is one of our most uncertain associations. 

(xxxv) ISOHDFS J223315-603224: This is the 



brightest source in our catalogue at 6.7^m and the 
second brightest at 15/im. It is a bright (K=13) M2V 
star, and is masked out of the GSFC CTIO catalogue, 
but its association with the EIS near-infrared cata- 
logue yields P^an < 0.001. 

In summary, we have found associations for 32 
out of the 35 ISO sources from Paper I. Of the remain- 
ing three, two (ISOHDFS J223256-603059 and ISO- 
HDFS J223256-603137) were detected only at 6.7pm 
and lie very close to sources which are very bright 
in that band, and which may well have compromised 
the background subtration procedure in that area, re- 
sulting in a false detection, or a significant shift in 
the source position. The third, ISOHDFS J223256- 
603513, was detected only at 15pm, and is a peak 
in an extended region of emission stretching from a 
bright source just below the southern boundary of 
our source detection region (see Figure 6 of Paper I) . 
Eight of the 32 identified sources are stars, leaving a 
total of 24 galaxies to be considered further. 



3 STAR FORMATION RATE 
ESTIMATORS 

In this Section we discuss and compare the various 
methods by which we might estimate the star forma- 
tion rates of our ISO sources. 

3.1 Methods for estimating star formation 
rates 

A multitude of methods have been advocated for the 
estimation of the star formation rates of galaxies, and, 
as discussed, for example, by Cram et al. (1998) and 
Granato et al. (2000), they agree to varying degrees. 
All attempt to detect the observational consequences 
of the formation of massive stars, and then use a 
model for the stellar initial mass function (IMF) to 
infer the total rate of formation of stellar mass. This 
inevitably introduces some uncertainty into the SFR 
estimates, due to the assumption of a universal IMF 
for all galaxies and the lack of a clear choice between 
competing models for it. More serious is the fact that 
the formation of massive stars in external galaxies, 
like that in our own Galaxy, is expected to take place 
in giant molecular clouds, enshrouded by dust, and, 
as argued, for example, by Silva et al. (1998) and 
Jimenez et al. (1999), the resultant emission from the 
star-forming region cannot be well represented by a 
simple model such as a dust-screen in front of a dust- 
free starburst. 

Despite these difficulties, there are several well- 
used prescriptions for estimating the SFRs of galaxies 
on the basis of simple broad-band fluxes, as these are, 
typically, all that are available for the large samples of 
galaxies used in censuses of the cosmic star formation 
history. The first of these is to use the U band magni- 
tude, arguing that the UV emission comes principally 
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Figure 5. Comparison of the results of applying to GRASIL model SEDs six star formation rate estimation methods 
from eqn. ^, namely those based on: (a) the 160nm luminosity, from Devriendt et al. (1999); (b) the 280nm luminosity, 
from Devriendt et al. (1999); (c) the U band luminosity, from Cram et al. (1998); (d) the 60^m luminosity from Cram 
et al. (1998) and Rowan— Robinson et al. (1999), denoted by the upper and lower sets of points, respectively; (e) the 3- 
1000//m integrated luminosity, from Devriendt et al. (1999); (f) the 1.4GHz luminosity, from Cram et al. (1998). For each 
method we plot the SFR estimate obtained by its application to the GRASIL SED models for MlOO (inverted triangle), 
M51 (diamond), NGC6946 (cross), NGC6090 (filled triangle), M82 (filled circle) and Arp220 (filled square) against the 
corresponding SFR value computed within the GRASIL model itself. 



from massive, young stars. There are two problems 
with this. Firstly, the U band (Acff = 365nm) is not 
far enough into the ultraviolet for the emission from 
older stars to be negligible, but it is possible to correct 
for this effect, by using model spectra to bootstrap 
from the U band to, say, 250nm, where the emission 
is dominated by massive stars. Cram et al. (1998) do 
this, obtaining the following relationship between the 
U band luminosity of a galaxy and its rate of massive 
(M > 5M0) star formation: 



SFR{M > 5Me 



L.(U) 



M, 



yr" 



1.5 X 1022 WHz~ 



(2) 



for an assumed ip{M) oc M~^'^ IMF running from 0.1 
to 100 Mq . An obvious objection to this is that it fails 
to account for the effects of dust in the star-forming 
galaxies, especially in the case that more massive star- 
bursts are dustier. The recent models of Devriendt, 
Guiderdoni & Sadat (1999) seek to model the spec- 
tral energy distributions (SEDs) of dusty starbursts 
from the UV to the submillimetre in a self-consistent 



fashion, taking account of the absorption of ultravi- 
olet Hght from massive young stars by the dust sur- 
rounding them, and its use in heating up the dust, 
leading to emission in the far-infrared and submil- 
limetre. Their models yield the following correlations 
between the star formation rate (for a Salpeter 1955 
IMF over [0.1,120] Mq) and the luminosities at two 
UV wavelengths: 



SFR 

Mq yr-i 

and 

SFR 

Mq yr-1 



ALA(280nm) 
7.7 X 1034W 



ALA(160nm) 
9.2 X 10^4^ 



(3) 



(4) 



Devriendt et al. (1999) caution against the over- 
interpretation of these best-fit correlations (about 
which their model galaxies scatter quite widely), but 
the non-linearity of the relationships they describe 
between SFR and UV luminosity is qualitatively what 
would be expected in a model in which the more mas- 
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sive a starburst is the dustier it is, and suggests that 
this effect may be corrected for, albeit crudely. 

The bolometric luminosities of luminous star- 
burst galaxies are dominated by thermal emission 
from dust in the far-infrared. In general, this emis- 
sion is a combination of extended ( "cirrus") emission 
from dust heated by the ambient interstellar radiation 
field, more localized "starburst" emission from the 
dust in regions of massive star formation heated by 
the UV flux from O and B stars, and, possibly, "AGN" 
emission from the hot dusty torus of an active galac- 
tic nucleus. The relative importance of these compo- 
nents varies between galaxies (e.g. Rowan-Robinson 
& Crawford 1989) , with the cirrus component declin- 
ing to higher luminosities, where ultra-luminous in- 
frared galaxies (ULIRGs) are seen to be mostly fu- 
eled by starburst emission, although some do have a 
dominant AGN component (Genzel et al. 1998): such 
simple modeling must be supplemented by considera- 
tion of emission from PAHs and inclusion of starburst 
components of different ages (Silva et al. 1998, Efs- 
tathiou et al. 2000) to fit the detailed SEDs revealed 
(e.g. Acosta-Pulido et al. 1996) by recent ISO spec- 
troscopy. Various attempts have been made to relate 
this far-infrared emission to the amount of ultraviolet 
light that must be absorbed by dust to produce it, by 
which route a galaxy's SFR (or, at least, the rate of 
formation of stars whose light is obscured by dust) can 
be estimated from its far-infrared luminosity. Rowan- 
Robinson et al. (1997) summarized previous work on 
this method in the relationship 



SFR 



M, 



1.5 X lO^^W 



(!) 



(5) 



where e ~ 1 is the fraction of the {vL,,) UV luminosity 
of the starburst that is re-emitted in the far-infrared, 
while (/) ~ 0(1) is a factor whose deviation from unity 
can account for variations of the IMF from the stan- 
dard Salpeter (1955) tpi^M) oc M''^-^^ form over the 
range 0.1 to 100 Mq : for example, it should be 3.3 
ff a Miller-Scalo IMF is preferred, and 1/3.1 ff the 
starburst forms only stars with M > I.GMq. Cram et 
al. (1998) presented a similar relationship (based on 
the ideas of Condon 1992), which reads 



SFR{M > 5Mg 



I/^(60/im) 



M, 



5.1 X 1023WHZ' 



(6) 



assuming the same ip^M) oc M~^'^ as in equation (j^. 

Similar expressions can be derived in terms of the 
integrated far-infrared/submillimetre emission. For 
example, the models of Devriendt et al. (1999) yield 
(for the same IMF as in equations ^ and ^ the rela- 
tionship 



SFR 



M, 



© yr 



Lir(3 - lOOOAtm) 



3.0 X 103«W 



(7) 



(cf. Eqns. ^ and ^ where the power law index deviates 
significantly from unity). 

As reviewed by Condon (1992) there is a well- 



known correlation between the far-infrared and deci- 
metric radio luminosities of actively star-forming 
galaxies, which is thought to arise from the fact that 
the bulk of the radio luminosity is produced by syn- 
chrotron emission from relativistic electrons spiraling 
in the remnants of supernovae originating in the same 
population of massive stars that produce the star- 
burst component to the far-infrared luminosity. Since 
the cirrus component to the far-infrared luminosity is 
not expected to have associated radio emission, it has 
been argued that the decimetric radio luminosity of a 
starburst galaxy may provide the cleanest handle on 
its star formation rate, and Cram et al. (1998) present 
the following form for that relationship (assuming the 
same IMF as before): 



L^(1.4GHz) 
4.0 X 1021 WHz" 



(8) 



SFRjM > 5Mq) 
Mq yr-1 

where, as before. Cram et al. consider only the for- 
mation of stars with 5 < M < IOOMq in an M"^-^ 
IMF. Additional SFR estimators exist, for example 
using the luminosity in the Ha line (e.g. Kennicutt 
1998), but these, too, are affected by dust, as shown 
by Rigopoulou et al. (2000) on the basis of near- 
infrared VLT-ISAAC spectroscopy of the Ha line in 
a subsample of the ISO-HDF-S objects considered 
here. 



3.2 Comparing star formation estimates 

It is natural to enquire how these SFR estimators 
compare, and this question was addressed empirically 
by Cram et al. (1998), using a somewhat heteroge- 
neous compilation of data (U band magnitudes. Ha, 
60/^m and 1.4 GHz radio fluxes) from a variety of 
sources. They found that the star formation rates de- 
duced from the integrated far-infrared and decimetric 
radio luminosities are well correlated over more than 
four orders of magnitude, but that the significant de- 
viations from linearity and greater scatter are seen for 
the relationships between the SFRs deduced from the 
1.4 GHz power and those from Ha and U band lu- 
minosities. Some observational effects (e.g. slit losses 
in Ha spectroscopy) may be contributing to these 
trends, but they are consistent with the qualitative 
expectations of the picture outlined above in which 
massive stars are formed in dusty environments, and 
reinforce the belief that this is the dominant mode of 
star formation in actively star-forming galaxies. 

Here we perform a complementary study, which 
attempts to circumvent the observational problems 
that inevitably affect the analysis of any heteroge- 
neous data sample drawn from many sources in the 
literature, by asking how these different star forma- 
tion estimators fare when applied to model SEDs for 
a range of types of star-forming galaxy. To do this we 
use the GRASIL models of Silva et al. (1998), which 
provide good fits to the UV-radio SEDs of six nearby 
galaxies, namely three starburst galaxies (Arp220, 
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M82 and NGC6090) with diflFering levels of activity, 
and three local spirals (MlOO, M51 and NGC6946): 
similar models of starbursts are presented by Efs- 
tathious, Rowan-Robinson & Siebenmorgen (2000), 
but they do not consider normal spirals, which is why 
we use the GRASIL models here. We refer the reader 
to the paper by Silva et al. (1998) for a discussion 
of the UV-mm SEDs, and to Silva (1999) for details 
of their extension into the radio, through consider- 
ation of separate thermal and non-thermal compo- 
nents to the radio emission, and note that Granato et 
al. (2000) have also compared star formation estima- 
tors through their application to GRASIL SEDs. 

Silva et al. (1998) quote the SFR value (averaged 
over the previous 5 x 10^ yr) corresponding to each 
SED model, for an assumed Salpeter (1955) IMF run- 
ning from 0.1 to 100 Mq. Hence, to compare the SFR 
estimators, by seeing how well each recovers those 
values, we must scale Eqns.(2) - (8) as required for 
them all to give the rate at which stellar mass would 
form given that reference IMF; this is analogous to the 
choice of the value of (f) in the formalism of Rowan- 
Robinson et al. (1997). The derivation of the correct 
scaling factor in each case requires consideration of 
the model through which the particular observed lu- 
minosity is related to the rate of formation of massive 
stars: in Eqns. (2), (3) & (4) the value of a particu- 
lar monochromatic UV luminosity is used as a direct 
tracer of these stars; underlying Eqns. (5), (6), & (7) 
is an assumption that some large (~ 1) fraction of 
their bolometric luminosity emerges [on a timescale of 
0{W^ yr)] in the far infrared (FIR), due to the re- 
processing by dust; and the method of Eqn. (8) relates 
the number of high-mass stars that become Type lb 
and Type II supernovae to the synchrotron emission 
produced by their remnants. 

In Appendix A, we show how each of these phys- 
ical models yields a scaling to be applied to Eqns. (2) 
- (8) to convert them to predictions for A4, the rate 
of formation of stellar mass, given our canonical IMF, 
which is a Salpeter (1955) IMF running from 0.1 to 
100 Mq, in our canonical cosmology, which is an Ein- 
stein - de Sitter universe with a Hubble constant of 50 
km s~^ Mpc~^. We find the following set of corrected 
versions of the estimators of eqns. (2) to (8): 

L,(f/)/2.6 X lO^^WHz-^ 
[AiA(280nm)/7.0 x lO^-'W] ^'^^ 
[AiA(160nm)/8.4 x lO^'-'W] ^'^"^ 
\Lx{60fJ.m)/1.5 X 10^*^ W 
L^(60Atm)/l.l X lO^^WHz-i 
[Lir(3 - 1000Aim)/2.3 x lO^^W] 
L^(1.4GHz)/6.9 x lO^^Wnz-V 



M 



Meyr- 



= < 



(9) 



Having made these various scalings we are ready 
to compare the SFR values obtained by applying the 
estimators of eqn. ^ to the six GRASIL SEDs with 
the values given for them by Silva et al. (1998), as 
shown in Figure ^ note that the results from the two 



60/im-based estimators are plotted together. The first 
thing to note from this figure is that the far infrared 
and radio estimators do far better at recovering the 
GRASIL SFRs than do the three UV-based ones. The 
integrated far infrared estimator of Devriendt et al. 
(1999) is the closest to the GRASIL SFR value for all 
models, while the 1.4 GHz is offset by a fairly con- 
stant factor of ~ 3, indicating that, while there may 
be a problem with its absolute normalisation, (pos- 
sibly caused by the way that the GRASIL SEDs are 
extended into the radio), this method works well for 
all the SED types considered, yielding accurate rela- 
tive SFR values for them. The 60^m-based estimators 
fare quite well, too, although the ratio of their SFR 
values to those from the GRASIL code are higher for 
starbursts (filled symbols) than for normal galaxies 
(empty symbols), suggesting that a far infrared colour 
term should be included in the relationship between 
SFR and far infrared luminosity. 

The three UV-based estimators display a much 
more complex behaviour across the range of galaxy 
types. The U band luminosity always gives a lower 
SFR estimate than the four integrated/far-infrared 
and radio estimators, and the discrepancy is appre- 
ciably larger for the starbursts than for the normal 
spirals, as expected if the UV light is generated in 
dustier environments in starbursts than in normal spi- 
rals. The two non-linear UV-SFR relations deduced 
by Devriendt et al. (1999) appear to over-estimate 
the star formation rates of the spirals, and one of 
the starbursts (NGC6090), but do reasonably well for 
M82 and, especially, Arp220. This could be because 
the Devriendt et al. (1999) models, although spanning 
the full range of galaxy types from inactive spirals to 
ULIRGs, are more directed at the understanding of 
starbursts and/or could simply refiect that the best- 
fit correlations yielding Eqns. (3) and (4) do not ex- 
press anything physical in the models. 

It is clear from these results that caution must be 
exercised in comparing SFRs in the literature, which 
might have been made using different estimators and 
with differing assumptions as to the stellar IMF and 
the exact specification of the astrophysical model un- 
derlying their application. Even when scaled to a com- 
mon reference IMF, the SFR estimators based on far 
infrared or radio luminosities are only consistent to 
within a factor of two, while the UV-based ones are 
seen to be far less secure. The consistency between 
the ratios of SFR values obtained by the four radio- 
and integrated/FIR-based estimators across the six 
SEDs is, of course, simply a manifestation of the uni- 
versality of the radio-FIR correlation (reviewed by 
Condon 1992), but our results do afi'ord some con- 
fidence that they are measuring quantities correlated 
with the "true" star formation rate, at least to the ex- 
tent that is well reproduced by the GRASIL models, 
which cannot be said for the UV-based estimators. 
On the basis of these results, the integrated IR lu- 
minosity appears to be the best SFR estimator out 



© 2002 RAS, MNRAS , 



Observations of HDF South with ISO - //. 13 



of the set we have investigated, so it is that which we 
shall use in what follows, when we deduce SFR values 
for our ISO sources. 



4 SPECTRAL ENERGY DISTRIBUTIONS 
AND STAR FORMATION RATES OF 
ISO-UDF-S SOURCES 

In this Section we estimate star formation rates for 
our ISO sources, through the application of the in- 
tegrated IR lu mino sity estimator judged to be the 
best in Section to GRASIL model SEDs fitted 
to the photometric data for each of the 24 galax- 
ies given in Table |l|, plus radio data (mostly upper 
limits) kindly provided in advance of publication by 
the ATNF HDF-S survey team (A. Hopkins, pnv. 
comm.). We do not expect our ISO sources to match 
exactly one of the six GRASIL models, but, to the ex- 
tent that the models span the range of SEDs of star- 
forming galaxies likely to feature in our mid-infrared 
survey, this method gives us a handle on the uncer- 
tainty in the deduced SFR value of each 750 source 
resulting from uncertainty in its true SED: note that 
the SEDs of the three normal spirals are very similar, 
and it is only when an appreciable starburst com- 
ponent kicks in that the shape of the SED begins 
to change significantly. To facilitate the choice of the 
GRASIL model most appropriate to each galaxy, we 
plot in Figure ^ the photometric data set for each 
source, together with the six model SEDs, normal- 
ized in each case to match the K band magnitude of 
the galaxy associated with the ISO-HDF-S source, 
or the I band magnitude in the one case (J223240- 
603141) of a galaxy lying outside the region of the EIS 
K band survey. From this figure we conclude the fol- 
lowing about the star formation rates of the 24 ISO- 
HDF-S galaxies: 

(i) ISOHDFS J223240-603141: The observed I 
band to 15/im colour suggests that this is a star- 
burst galaxy, rather than a normal spiral dominated 
by cirrus emission: the lack of a detection at 7/^m is 
just consistent with that, and J223240-603141 does 
lie close to a very bright 7/im source, so it is possi- 
ble its 7/im flux estimate has been corrupted by the 
subtraction of the negative lobes of the bright source. 
The lack of detections at 1.4, 2.5 and 4.9 GHz ar- 
gue against the starburst being as extreme as that 
in Arp220, and the optical colours favour an SED like 
that of NGC6090 over that of M82. Using the 1.4GHz, 
60/im (Cram et al.), 60/im (Rowan-Robinson et al., 
1997, with e = (j) = 1), 3-1000/im prescriptions we 
obtain SFR estimates of 43, 34, 12 and 20 M© yr~\ 
respectively, from adopting this SED. On the basis of 
Section pj, we adopt the last of these, so our best es- 
timate of the true SFR for this galaxy is 20Mq yr~^. 

(ii) ISOHDFS J223243-603242: This is a 
slightly ambiguous case. The optical/near-infrared 
colours of this galaxy fit the M82 model SED very 



well, as does the 7^m flux, once the model has been 
normalized to the observed K band flux of J223243- 
603242; however the 15/im flux seems a little low and 
the lack of a 1.4 GHz detection is marginally incon- 
sistent with the M82 model. Using the M82 model 
SED, the same four SFR estimators give 145, 121, 
43, 60 Mq yr"\ while, if we had adopted an MlOO 
SED instead, we would have obtained values of 22, 
9, 3 and 7Mq yr"\ respectively. The M82 SED is a 
better fit overall, so we adopt a best guess SFR value 
for this galaxy of 60 Mq yr~^, noting that, while the 
precise SFR value is uncertain, it is clear this galaxy 
is forming stars at a rate of several tens of Mq yr~^. 

(iii) ISOHDFS J223243-603351: This is one of 
the four ISO-HDF-S sources detected in the radio 
and the strength of these detections (at 4.9 and 8.6 
GHz), together with the optical/near-infrared SED 
of the galaxy, which rises into the UV, suggests that 
this may be an AGN, and this is confirmed by the 
presence of a broad line in its optical spectrum: we 
therefore, do not estimate the SFR in this source. 

(iv) ISOHDFS J223243-603441: This is a sec- 
ond source, like J223243-603242, which has an 
optical/near-infrared SED well matching the star- 
burst models, but ISO fluxes which do not unam- 
biguously support that interpretation, as the 15pim 
flux is lower than would be expected for a starburst 
on the basis of its SED up to 6.7/im. From Fig. ^ we 
see that this source is close to a much brighter source 
at 15/im, so it is possible that its flux in that band 
has been under-estimated, due to the difficulty of sub- 
tracting the negative lobe from the brighter source. If 
we assume that, and adopt an NGC6090 SED, then 
we obtain SFR values of 216, 169, 61, 108 Mq yr"^ 
from the usual four estimators, while, if we take the 
6.7/xm-15/xm colour at face value, favouring an M51 
SED, we obtain, instead, values of 49, 28, 10 and 19 
Mq yr~^. As a compromise, we adopt a best guess of 
50 Mq yr~^, but note that is is uncertain by a factor 
of two, at least. 

(v) ISOHDFS J223244-603110: The optical to 
near-infrared SED of this galaxy is a little unusual, 
suggesting a mis-match in the apertures used to mea- 
sure the optical and near-infrared magnitudes. Fixing 
the model SEDs to match the K band magnitude we 
obtain a good fit to the M82 and NGC6090 models, 
yielding an SFR estimate of 4 Mq yr~^ from applying 
the integrated IR luminosity estimator to the fitted 
M82 SED mode, so this is not a very strong starburst. 

(vi) ISOHDFS J223244-603455: The proximity 
of a brighter 7pm source might explain the lack of a 
detection in that band, due to the lobe-subtraction 
problem, but it seems more likely, on all evidence, 
that this is a normal spiral and not a starburst. Using 
the NGC6946 model SED the integrated Lir estima- 
tor yields 4Mq yr~^, a modest star formation rate, 
consistent with this interpretation. 

(vii) ISOHDFS J223245-603226: Again, the 
lack of a 7/im detection slightly confuses an other- 
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Figure 6. Spectral energy distributions for the 24 objects from Table |l| identified as galaxies. The squares mark the 
U V/optical/near-IR and radio photometric data for the objects, from Table ^ while the solid, dotted, dashed, dash- 
dot, dash-dot-dot-dot and long— dashed lines show, respectively, the GRASIL model fits to the SEDs of Arp220, MlOO, 
M51, M82, MGC6090 and NGC6946, redshifted as appropriate for each galaxy The ISO data from Table 8 of Paper 
I arc plotted as error bars or upper limits, as appropriate. High resolution version of this figure available from 
astro.ic.ac.uk/hdfs. 



Figure 6 - continued Spectral energy distributions for the 24 objects from Table ^ identified as galaxies. The squares 
mark the UV/optical/near-IR and radio photometric data for the objects, from Tab e ^ while the solid, dotted, dashed, 
dash-dot, dash-dot-dot-dot and long-dashed lines show, respectively, the GRASIL model fits to the SEDs of Arp220, 
MlOO, M51, M82, MGC6090 and NGC6946, redshifted as appropriate for each galaxy The ISO data from Table 8 of Paper 
I are plotted as error bars or upper limits, as appropriate. High resolution version of this figure available from 
astro.ic.ac.uk/hdfs. 



wise fairly confident identification of this witfi an 
M82- or NGC6090-type starburst, and, once more, 
this source is close to a strong Jfim source, so this may 
be due to the lobe-subtraction problem. If we adopt 
the NGC6090 SED model, which gives a slightly bet- 
ter fit to the optical/near-infrared SED, we obtain a 
best guess SFR of SGAf© yr~^ for this galaxy. 

(vui) ISOHDFS J223245-603418: This is one of 
the most unambiguous starburst detections. Its UV- 
15/xm colours are well fitted by the M82 SED (and 
almost as well by the Arp220 model), but the solid 
radio detections at 1.4, 2.5 and 4.9 GHz suggest a star- 
burst stronger than M82. Using the M82 and Arp220 
SEDs we obtain SFRs of 40 and SOOM© yr"\ using 
the integrated Lir, estimator, confirming that this is 
a powerful starburst, although leaving some doubt as 
to its true star formation rate: we adopt a best guess 
of 100 Mq yr~^, but note this is uncertain by a factor 
of two, at least. 

(ix) ISOHDFS J223247-603335: Another case 
where an unambiguous discrimination between spi- 
ral and starburst is difficult: the ISO fluxes are a lit- 
tle lower than expected for the starburst fit to the 
optical/near-infrared SED. Taking the NGC6090 and 
MlOO SEDs, which bracket the data points, we obtain 
star formation rates of 130 and 16 Mq yr~^, from ap- 
plication of the integrated Lir estimator, so there is 
quite some uncertainty in the SFR for the source: we 
adopt a value of 50Mq yr~^, again with a factor of 
two uncertainty, at least, although it is clear that this 
galaxy is definitely forming stars actively. 

(x) ISOHDFS J223251-603335: The SED of 
this source is quite well fit by the NGC6090 model, 
with which the integrated Lir estimator yields an 
SFR estimates of 10 Mq yr~^, so it is clear this is 
a modest starburst. 

(xi) ISOHDFS J223252-603327: This radio 
source, detected solidly at 1.4 GHz, is at the centre of 
extended emission in both ISO bands, which may ex- 
plain why the mid-infrared fiuxes slightly exceed the 
predictions of the starburst models. The radio detec- 
tion is in excellent agreement with the Arp220 model, 
and if we adopt that we obtain 500 Mq yr~^, using 
the integrated Lir, estimator. It was noted in Subsec- 



tion 2.3 that this source becomes more point-like at 
longer wavelengths, as one passes through the optical 
into the near-infrared, so it could contain an obscured 
AGN, and, hence, the star formation rate estimated 
here may be an over-estimate. 

(xu) ISOHDFS J223254-603115: As noted in 
Section 2.3, there appears to be a mismatch in the 



apertures used for the optical and near-infrared pho- 
tometry of this galaxy. If we normalize the GRASIL 
models to the K band magnitude of this galaxy, then 
we see that its SED is in good agreement with the 
NGC6090 and M82 models, with the Arp220 model 
marginally disfavoured by the lack of a radio detec- 
tion. Applying the integrated Lir estimator to the 
NGC6090 model we obtain an SFR of IIM© yr-\ 
which we adopt as our best guess; if, instead, the 
M82 model were used, that figure would be negligibly 
higher, at 12 Mq yr~^. 

(xni) ISOHDFS J223254-603129: From Fig- 
ure ^ we see that the ISO fiuxes of this source are in 
good agreement with the three starburst model SEDs, 
when they are normalized using the galaxy's K band 
magnitude. If we apply the integrated Lir estimator 
to the NGC6090 model, which gives the best over- 
all fit, we obtain a best guess SFR estimate of 1.2 
Mq yr-\ 

(xiv) ISOHDFS J223257-603305: The 15^m to 
near-IR colour of this galaxy implies that it is a star- 
burst, and the lack of radio detections excludes a 
burst as extreme as Arp220. Taking the NGC6090 
model, we obtain an SFR estimate of 43 Mq yr~^, 
using the integrated Lir estimator, and we adopt this 
as our best guess. 

(xv) ISOHDFS J223302-603323: Figure | 
shows that this source is securely identified as a cirrus 
galaxy, rather than starburst, and, adopting the MlOO 
model, we obtain a best guess SFR of 10 Mq yr~^, 
using the integrated Lir estimator. 

(xvi) ISOHDFS J223303-603336: The 7^m to 
near-IR colour of this source suggests that it is a star- 
burst, rather than a cirrus galaxy, but the lack of 
a 15/im detection is puzzling, in that case, although 
once again it may be due to uncertainty in the sub- 
traction of the negative lobes from a nearby source 
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Figure 6 - continued Spectral energy distributions for tlie 24 objects from Table gj identified as galaxies. The squares 
mark the U V/optical/near— IR and radio photometric data for the objects, from Table ^ while the solid, dotted, dashed, 
dash— dot, dash— dot— dot— dot and long-dashed lines show, respectively, the GRASIL model fits to the SEDs of Arp220, 
MlOO, M51, M82, MGC6090 and NGC6946, redshifted as appropriate for each galaxy. The ISO data from Table 8 of Paper 
I are plotted as error bars or upper limits, as appropriate. High resolution version of this figure available from 
astro.ic.ac.uk/hdfs. 
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that is bright in both 750 bands. The Arp220 SED 
is ruled out by the lack of radio detections, and, if 
we adopt the NGC6090 model, we obtain an SFR 
of 16Mq yr~^, using the integrated Lir estimator, 
while, if the MlOO model is preferred, this falls to 
2Mq yr~^. We shall adopt a best guess of 5 Mq yr~^, 
but note that this is one of our most uncertain SFR 
estimates. 

(xvu) ISOHDFS J223306-603349: From Fig- 
ure ^ we see that this galaxy was detected in the ra- 
dio at 1.4 and 2.5 GHz at almost exactly the level ex- 
pected for the GRASIL NGC6090 SED normalized to 
the galaxy's observed K band magnitude, but that the 
two ISO fluxes look a little low to match that inter- 
pretation, especially the 15/xm flux. The radio detec- 
tions are a strong indication that this is a starburst, 
rather than a cirrus galaxy, however, so we do adopt 
the NGC6090 SED, which yields a best guess SFR 
of 60Mq yr~^, from the integrated Lir estimator, al- 
though we record that, if we had adopted a NGC6946 
SED, this would have fallen to 10 AIq yr^^. 

(xviii) ISOHDFS J223306-603436: The mid- 
IR colour of this source suggests that it is a starburst, 
rather than a cirrus galaxy, and its optical/near-IR 
colours flt an (appropriately normalized) M82 model 
all the way from the B to K band. Adopting that SED 
yields a best guess SFR of 74±15 Af© yr"^ 

(xix) ISOHDFS J223306-603450: As illus- 
trated in Figure ^ the 15/im/near-IR colour of this 
source is that of a starburst, not a cirrus galaxy, but 
the extant data are not able to distinguish between 
an Arp220-type galaxy and a more modest burst, like 
M82, with the lack of detections at 7fim and in the 
radio consistent with both possibilities. The Arp220 
SED would yield an SFR of 130 Mq yr"\ while the 
M82 model gives 17Mq yr~^, so we shall adopt an ad- 
mittedly uncertain compromise flgure of 50i\f© yr~^ 
as our best guess SFR. 

(xx) ISOHDFS J223307-603248: We cannot 
unambiguously distinguish whether this source is 
a starburst or cirrus galaxy. The 7/15/im colours 
favours the former, as does the shape of its 
UV/optical/near-IR SED, and, if we adopt the 
NGC6090 SED we obtain an SFR estimate of 



65Mq yr~ , while a lower value of 10 Mq yr~ would 
result from using an MlOO model instead: we shall 
adopt the former value as our best guess, but note 
that we cannot confldently exclude an SFR value a 
factor of four lower. 

(xxi) ISOHDFS J223308-603314: This source 
is solidly identified as a starburst, as its SED fits the 
GRASIL NGC6090 model from the U band all the 
way to 15/xm. With that SED we estimate its SFR to 
be 46Mq yr~^, using the integrated Lir estimator, 
(xxu) _JSOHDFS J223312-603350: As noted in 
the redshift of this galaxy is highly un- 
hence, so is its star formation rate. The 



Section 2.3, 
certain, and, 

shape of its UV/optical/near-IR SED suggests that it 
is a starburst, as, more strongly, does its mid/near-IR 
colour, but it is difficult to judge the strength of its 
burst from the extant data. We adopt an NGC6090 
SED to obtain a best guess SFR estimate of 15 
Mq yr~^, although we note that choosing an Arp220 
model instead would have yielded 130 Mq yr~^. 

(xxui) ISOHDFS J223312-603416: This is an- 
other source with a very uncertain redshift, based only 
on JHK photometry, and, furthermore, its 7fim flux is 
highly uncertain, due to its being located in a region 
of extended 7fj,m emission, as shown in Figure ^ The 
15/^m/near-IR colour of this source suggests that it is 
a starburst, and the Arp220 and NGC6090 SEDs give 
SFR estimates of 370 and 46 Mq yr~^, respectively. 
We adopt a value of 100 Mq yr~^ as our best guess, 
but noting that, in addition to the uncertainty due to 
the poorly constrained SED, there is also uncertainty 
over the redshift of this galaxy. 

(xxiv) ISOHDFS J223314-603203: Another 
source whose mid/near-IR colour indicates it to be 
a starburst, rather than a cirrus galaxy, but for which 
we are unable to assess the strength of the burst given 
the current data. The Arp220 and NGC6090 SED 
models give SFR estimates of 9 and 1 Mq yr~^, re- 
spectively, indicating either way that this is a modest 
starburst: we shall adopt a value of 3 Mq yr~^ as our 
best guess, noting that Section 3 showed this to be 
one of our most uncertain associations. 
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Figure 6 - continued Spectral energy distributions for the 24 objects from Table gj identified as galaxies. The squares 
mark the UV/optical/near-IR and radio photometric data for the objects, from Table ^ while the solid, dotted, dashed, 
dash-dot, dash-dot-dot-dot and long-dashed lines show, respectively, the GRASIL model fits to the SEDs of Arp220, 
MlOO, M51, M82, MGC6090 and NGC6946, redshifted as appropriate for each galaxy. The ISO data from Table 8 of Paper 
I arc plotted as error bars or upper limits, as appropriate. High resolution version of this figure available from 
astro.ic.ac.uk/hdfs. 
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Figure 7. The distribution of star formation rates for 
the 22 ISO HDF— S sources we believe to be star form- 
ing galaxies, omitting the galaxy with the highest SFR 
(ISOHDFS J223252-603327 at SOOM© yr"!), which may 
be harbouring an obscured AGN. 




Figure 8. A comparison of the cumulative redshift dis- 
tributions of /SO-HDF-S galaxy sample (solid line) with 
those from three HDF-S photometric redshift catalogues; 
the dot-dot-dot-dashed line is that of MRR, the dashed 
line that of Gwyn, and the dotted line that of the SUNY 
group. 



5 THE CONTRIBUTION OF 

/50-SELECTED SOURCES TO THE 
STAR FORMATION HISTORY OF THE 
HDF-S 

In Figure ^ we plot the distribution of SFR values 
for 22 /SO-HDF-S sources we believe to be star- 
forming galaxies, omitting the galaxy with the high- 
est SFR (ISOHDFS J223252-603327 at 500Mq yr"^), 
which may be harbouring an obscured AGN. As can 
be seen from this, the ISO-HDF-S galaxies are typi- 
cally active starbursts (median SFR — 43 Mq yr~^), 
although they do cover a range of two orders of 
magnitude in star formation rate. These values are 
computed under the assumption of an Einstein - 
de Sitter Universe with a Hubble constant of Ho — 
50 km Mpc~^. Using the slightly non-linear in- 
tegratd Lir estimator of Devriendt et al. (1999) we 
have the SFR varying with luminosity distance, ^l, 
as di,'^! and, hence, with Hubble constant as _ffo'**. If, 
instead, we had assumed a cosmology with 51m =0.7 
and JIa = 0.3, we would have deduced, for the same 
value of _ffo SFR values that were factors of (11, 20, 
47 & 78) per cent higher for galaxies at redshifts of 
(0.1, 0.2, 0.5 & 1.0) respectively. 

For five of our galaxies, SFR estimates have been 
made by Rigopoulou et al. (2000), on the basis of 
Ha luminosity and also using a far-infrared (FIR) 
estimator (Franceschini et al., in preparation) which 
makes use of the 15jUm flux and assumed a far- 
infrared/mid-infrared luminosity ratio of ~ 10 as ap- 
propriate for a galaxy with a spectral energy distri- 
bution like that of M82. Rigopoulou et al. (2000) as- 
sumed a cosmology with flu = 0.7 and SI a = 0.3 
and a Salpeter (1955) IMF over the mass range 
[1,1OO]M0, so, before we can compare their results 
with ours, we must scale their SFR values to account 
for these effects, using methods discussed in Appendix 
A. Under our canonical assumptions, their FIR-based 
estimator gives SFR values of (113, 89, 69, 220, 
65) Moyr"^ for galaxies (J223245-603418, J223245- 
603226, J223247-603335, J223252-603327, J223257- 
603305), respectively, for which our adopted values 
are (100, 36, 50, 500, 43)Af0yr~^ , which is a rea- 
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sonable level of agreement, given the uncertainties in 
the true SEDs of these sources. With the appropri- 
ate scalings, the Ha-based estimator of Rigopoulou 
et al. (2000) gives SFR raw values of (2.1, 6.0, 5.9, 
74, 12) Mgyr"^ for these five galaxies, which are fac- 
tors of 5-50 smaller than those estimated from the 
FIR. Rigopoulou et al. (2000) argue that, on the ba- 
sis of the V-K colours of these galaxies, one would 
deduce an extinction correction of only ~ 4 to the 
Ha SFR, suggesting that optical data alone are not 
sufficient to derive a good SFR value, even when they 
do lead an extinction correction; this is in accordance 
with the work of Silva et al. (1998) and Jimenez et 
al. (1999), mentioned above, who argue on theoretical 
grounds that the emission from star-forming regions 
cannot be well represented by a simple model of a 
dust-screen in front of a dust-free starburst. 

5.1 Computing the raw star formation rate 
density 

We may use these SFR estimates to derive constraints 
on the star formation history of the Universe, at least 
in the redshift interval in which our sources are found, 
i.e. z < 0.6 for the most part. In fact, as Figure H 
shows, the redshift distribution of the galaxies asso- 
ciated with our ISO sources is consistent with that of 
similarly bright (i.e. / < 22.5) galaxies in the HDF- 
S as a whole, at least as judged from photometric 
redshift catalogues produced for the field: two-sided 
Kolmogorov-Smyrnov tests reveal that the cumula- 
tive redshift distribution of the ISO IDs has proba- 
bilities of 0.30, 0.43 and 0.75 of being drawn from the 
same population as that yielding the n{z) distribu- 
tions found in the photometric redshift catalogues of 
MRR, Gwyn and the SUNY group, respectively. 

We divide this region into two equal volume bins 
{0 < z < 0.43 and 0.43 < z < 0.6), in which we com- 
pute the star formation rate density, via the available 
volume technique. The volume in which each galaxy 
could have been observed within any given redshift 
slice was estimated by considering its luminosity and 
best fit SED and hence determining the effective area 
over which it could have been observed at given red- 
shift using the area as a function of 15/im fiux lim- 
its (Figure 13 of Paper I). The contribution of each 
galaxy to the global star formation rate density, p», 
was then calculated as the ratio between the star for- 
mation rate determined for that galaxy and its avail- 
able volume, and the sum of these individual contri- 
butions was then taken over all galaxies in each bin. 
Note that there exists one galaxy (ISOHDF J223303- 
603336, which may be affected by the lobe subtraction 
problem) lying the range z < 0.6 which does not have 
a detection at 15/xm. Whether this galaxy is excluded, 
or included, using an available volume computed with 
an estimate of what its 15/im flux "should" be, on the 
basis of the best fitting GRASIL SED, makes a negli- 
gible difference to the raw value of p, obtained, so, for 



definiteness, we neglect it, leaving us with a 15^m- 
selected sample. 

This yielded values of p, — 0.03 ± 0.2dex (i.e. 
p. = 0.03to;oi) and p, = 0.07 ± O.ldex {i.e. p, = 
0.07tom) Meyr-^Mpc"^ for < 2 < 0.43 and 
0.43 < 2 < 0.6, respectively, where we have assumed 
an Einstein - de Sitter Universe, star formation ac- 
cording to a Salpeter (1955) IMF over the mass range 
[0.1,100] Mq and taken our SFR estimates for indi- 
vidual galaxies to be uncertain by a factor of two. In 
Figure y we show a comparison of these results and a 
compilation (principally that of Haarsma et al. 2000) 
of similar constraints derived recently from a variety 
of methods in a number of wavebands, all using the 
same assumed cosmology and IMF. One notable fea- 
ture of this is that our raw results are in excellent 
agreement with another recent study using ISO data, 
that of Flores et al. (1999) based on their study of the 
CFRS 14154-52 field. 

5.2 Sampling variance in our p, results 

The relatively small volume of our survey region 
means that our star formation rate density values 
are subject to a relatively large sampling variance. 
The volume of a cone < z < 0.43 with solid an- 
gle 20 sq. arcmin (equal to our survey area) is ~520 
h~^Mpc'^, while a frustum with the same solid an- 
gle and 0.43 < 2 < 0.6 is ~570 /i"^Mpc^. We can 
thus consider each survey section to have a volume 
of ~ 550/i~'^Mpc^; cubic cells of the same volume 
would have sides of 8.2 Mpc, while spheres would 
have radius of 5.0 Mpc. These volumes are small 
enough that the variation in mean density between 
surveys carried out in different regions of the Universe 
are expected to be significant. Using the power spec- 
trum of Peacock & Dodds (1994) and the conversion 
from length scale to effective wavenumber for cubic 
cells (Peacock 1991) we estimate that the fractional 
density fluctuations in cubic cells of this volume at 
redshift zero will be = 1.05, while that of spheres 
would have cr^ — 1.15. Note that the variances in local 
cells with the same geometry as our survey volumes 
are likely to be larger than this, as their elongated 
shape would lead to a window function that lets in 
more small-scale power than spherical or cubic win- 
dow functions do. Several more effects further com- 
plicate an assessment of the sampling variance in our 
SFR density estimates. Firstly, there is evolution in 
the density field with redshift across the two sam- 
pling volumes. At z = 0.2 and z — 0.5, in the middle 
of our two redshift bins, the linear theory variances 
in a given volume would be factors of (1 -I- z)'^, i.e. 1.4 
and 2.3, respectively, below their values at redshift 
zero, while non-linear effects as cr —> 1 are likely to 
make the decrease more pronounced than these lin- 
ear theory estimates suggest. Furthermore, the shape 
of our survey volumes means that they have more 
volume at higher redshift, where the density field is 
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Figure 9. A compilation of constraints on the star formation history of the Universe, mainly taken from Haarsma et al. 
(2000): in all cases, an Einstein — de Sitter universe, with a Hubble constant of 50 km s~^ Mpc~^ was assumed, as was star 
formation taking place with a Salpeter (1955) IMF, over the mass range [0.1, 1OO]M0, and the dust extinction corrections 
of Steidel et al. (1999) were applied. The symbols are as follows: asterisk - Hughes et al. (1998); filled upward-pointing 
triangle - Flores et al. (1999); empty square — Connolly et al. (1997); plus sign - Gallego et al. (1995); empty downward- 
pointing triangle — Steidel et al. (1999); five-pointed star — Madau et al. (1996); diamond - Treyer et al. (1998); filled 
downward-pointing triangle - Rowan-Robinson et al. (1997); empty upward-pointing triangle — Lilly et al. (1996); empty 
circle — Haarsma et al. (2000); filled square - this work, with confidence intervals computed according to the sampling 
variance treatment of Appendix B, assuming cr^ = 1. 



less evolved. A further complicating factor is the bias 
between the galaxy and mass distributions, which is 
likely to lead to a higher sampling variance for the 
total SFRs in our two redshift bins than that for the 
meiss they contain, and an additional uncertainty in 
the redshift variation across the bins: KaufTmann et 
al. (1999) show that the evolution of galaxy clustering 
strength over < z < 1 is a function of cosmology, at 
least in their galaxy formation models. Clearly, the 
quantitative assessment of these various factors is be- 
yond the scope of the current paper, so we shall base 
this study of the effects of sampling variance on p» 
estimates on the assumption that the variance in the 
mass contained in our two redshift bins is cr'^ ~ 1 and 
that = 0.5 and = 2 are likely to be conservative 
bounds to the true values. 

The variance, is equal to o'^/p'^, where is 
the variance in the mass density field, and p is its 
mean value. A crude approach would be then to argue 
that 

cr = — ~ = (Jlnp, (10) 

P P 

in which case (Jiog^gp = log^o(e)o-. This would then 
mean that the estimates for in our lower (higher) 



redshift bins would have sampling errors of ±0.52dex 
(±0.34dex), respectively, so that the sampling vari- 
ance in p* exceeds the uncertainty from the individual 
SFR estimates for both bins. 

This method neglects the fact that, since ~ 1, 
the probability distribution function (PDF) for the 
cosmological density field will have been significantly 
skewed by gravitational evolution. This skewness will 
be refiected in strongly asymmetric error bars on our 
p* estimates due to sampling variance: with such a 
skewed PDF, our randomly-selected survey volume 
(the selection of the HDF-S region was based on the 
presence of a z ~ 2 quasar, which should have no 
bearing on the properties of galaxies at z < 1) is 
much more likely to be underdense than overdense, 
and, hence, we are much more likely to have mea- 
sured a value of p* that is lower than the global mean 
than one higher that it. Estimating the level of such 
an effect is very difficult. A full analysis would require 
a large numerical simulation of the cosmological den- 
sity field, coupled to a galaxy formation model capa- 
ble of predicting accurately the sites of formation of 
the class of galaxies detected in our ISO survey, which 
is far beyond the scope of this paper. Simple analytic 
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models prescriptions do exist for following the grav- 
itational evolution of the PDF, but most are based 
on approximations {e.g. perturbation theory or the 
Zel'dovich approximation) that break down by a ~ 1. 
Bernardeau & Kofman (1995) have, however, shown 
that a lognormal model continues to give a good fit to 
the PDF derived from Af-body simulations of a Cold 
Dark Matter universe until at least a 1.5, so this is 
one analytic form that might be used. 

In Appendix B we show how this approxima- 
tion may be used, within a Bayesian framework, to 
compute the 68 per cent confidence intervals for , 
the global star formation rate density at the redshifts 
corresponding to the centres of our two bins, given 
our raw p, values and our estimate that the mat- 
ter variance in cells of size equal to our survey vol- 
umes is 0.5 < cr'^ < 2. This analysis reveals that, for 
a best guess of = 1, 0.62 < p»/ p* < 3.2, while 
increasing or decreasing by a factor of 2 barely 
changes the lower limit to Pt/pt, while the upper 
limit increases to 5.0 or decreases to 2.3, respec- 
tively. Thus we see that sampling variance introduces 
a larger uncertainty into the estimation of p, than 
that caused by the uncertainties in the SFRs of the 
individual galaxies in our survey. Assuming — 1, 
we estimate 0.02 <p^< 0.10 and 0.04 <}, < 0.23 
Mgyr-^Mpc"^ for the z < 0.43 and 0.43 < z < 0.6 
bins, respectively, where we define confidence inter- 
vals solely using our sampling variance analysis. 

We plot these confidence intervals on Fig. ^. 
These results are model-dependent, and neither the 
simple variance estimation nor the lognormal model 
of Appendix B is satisfactory, but both indicate that 
there is a large uncertainty associated with p, esti- 
mates determined from volumes as small as our ISO 
survey of the HDF-S. Our results are consistent with 
those of the /50-based results of Flores et al. (1999) 
(whose CFRS field is an order of magnitude larger, 
and so should yield p, values with much lower sam- 
pling variances) and with the radio-based results of 
Haarsma et al. (2000), but it is clear that, given the 
large sampling variance inevitable for such a small 
survey volume, they cannot place tight constraints on 
the star formation history of the Universe by them- 
selves. 



6 DISCUSSION AND CONCLUSIONS 

In this paper and its companion (Paper I) we have 
presented results from our ISO survey of the Hubble 
Deep Field South. Here we sought optical identifica- 
tions for the ISO sources found in Paper I, obtain- 
ing reliable associations for 32 out of the 35 sources: 
these associations should be more secure than those 
made by Mann et al. (1997) using our initial analy- 
sis of our corresponding ISO survey of the northern 
Hubble Deep Field, thanks to the inclusion of cor- 
rections for CAM image distortions, which were not 



well characterised in 1997. Of these 32 sources, a to- 
tal of twenty two were identified as spiral or star- 
burst galaxies, eight of which have spectroscopic red- 
shifts (from the work of Rigpoulou et al. 2000 and 
Glazebrook et al. 2002) and the remaining fourteen 
of which have had photometric redshifts estimated 
by ourselves and others. We estimate that our pho- 
tometric redshifts should be accurate to 5z = 0.1 or 
so, on the basis of the set of eight objects for which 
spectroscopic redshifts are known, although we have 
noted individual cases where the error is likely to be 
greater than that, for example when a redshift has 
been estimated solely on the basis of JHK photom- 
etry. We found that the redshift distribution of the 
galaxies associated with our ISO sources is consistent 
with that for similarly-bright optical galaxies in the 
HDF-S region as a whole. 



We reviewed a series of methods commonly used 
to determine star formation rates for galaxies in sur- 
vey data, typically comprising a single broad-band 
flux for each object. We assessed the ability of these 
to reproduce the SFRs of models of actively star- 
forming galaxies, finding, as others (e.g. Granato et 
al. 2000) have, that those probing their far-infrared 
emission fared much better than those based on UV 
luminosity, indicating that massive star formation in 
these galaxies takes place in dusty regions: the work of 
Rigopoulou et al. (2000) shows further that this dust 
obscuration is not well corrected for using optically- 
determined extinction values. All these methods look 
for indications of the formation of massive stars, so a 
further complication in the determination of the ab- 
solute rate at which stellar mass is being formed in 
a given galaxy is uncertainty in the IMF, since the 
vast majority of the mass resides in stars not directly 
probed by these methods, so SFR estimates differing 
by factors of a few can result if differing IMFs are 
assumed. Caution must be exercised when applying 
these estimators in situations where the SED type of 
the galaxies under study are constrained only by a 
small number of broad-band fiuxes, suggesting that 
constraints on the cosmic star formation history are 
best obtained in well-studied fields, such as the Hub- 
ble Deep Fields, where rich, multiwavength datasets 
are available. The small areas of the Hubble Deep 
Fields do, however, mean that significant sampling 
variances exist for estimates of p* at 2 < 1, and we 
showed, via two simple methods for assessing their 
magnitude, that while our results are consistent with 
those of previous authors, notably Flores et al. (1999) 
and Haarsma et al. (2000), they cannot by them- 
selves yield tight constraints on the star formation 
history of the Universe, due to these sampling ef- 
fects. Further details of this project can be found at 
astro . ic . ac .uk/hdf s. 
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APPENDIX A: CONVERSION OF SFR 
ESTIMATORS TO A CANONICAL IMF 

As mentioned in Section 3.2, there are three types 
of physical model underlying the SFR estimators of 
Eqns. (2) - (8). It follows, by the definition of M, the 
rate of formation of stellar mass (in M0yr~^), that 



M 



MiPiMJdM, 
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Ml 



where tp{M) is the IMF , and stars are being formed 
over the mass range [Ml, Mu]- 

The first class are those (from Eqns. 5, 6 & 7) 
based on the assumption (Thronson & Telesco 1986, 
Rowan-Robinson et al. 1997) that some large (e ~ 1) 
fraction of the bolometric luminosity, Lboi, generated 
in any burst of star formation emerges in the far in- 
frared (FIR), due to the reprocessing by dust of the 
light from young stars: i.e. 
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where L{M) is the luminosity of a newly-formed star 
of mass M (and it is assumed that only main sequence 
stars contribute) and that such stars contribute to 
Lfih for time tFiR(M). The crucial timescale here is 
the time taken for the dust cloud around the young 
star to be destroyed and that is assumed (Thronson 
& Telesco 1986) to be independent of mass, and to be 
TpiR ~ 0(10^ yr), in which case 



Mu 

/ L{M)ip{M)AM, 

'Ml 



(A3) 



Al 



and AS 



so that combining Eqns 
M oc Lfir ■ 



Q L{M)^{M)dM 

M 
'7 ' 



(A4) 



where M and Zboi are, respectively, the average mass 
and bolometric luminosity of stars formed according 
to the particular IMF: i.e. 



M : 



and 



Im" Mi>{M)dM 
CHM)dM 



/"^^ LiM)i,{M)dM 
k"V'(M)dM 



(A5) 



(A6) 



So, we have shown, following Thronson & Telesco 
(1986) that the SFR (in Moyr'^) per unit of far in- 
frared luminosity (say, at 60/im, or integrated over 
3-1000/im) is proportional to A//Lboi. 

Several choices exist for the form of L{M) to use. 
Telesco & Gatley (1984) assume a double power law 
form, namely L{M)/Lq = A ■ {M/Mq)°', where 



{A,a) = 



(1.3,3.6) 0.1 < M/Mq < 10 
(8.1,2.8) 10 < M/Mq < 60, 



(A7) 



using which yields the following set of M/ Lboi values 
for the tp{M) oc M~^ models of interest here: 



M/L 



bol 



W-nie/Le 
fl.33 : {Ml,Mu) 



(O.l,lOO)M0,a; = 2.35 
h.m : {Ml,Mu) = (0.1,120)Afo,2: = 2.35(A8) 
[2.52 : (Ml,Mu) = (O.l,lOO)A/0,a; = 2.50, 

reproducing the results in Table 3 of Thronson & Tele- 
sco (1986), from which Rowan-Robinson et al. (1997) 
derived the value that (f} should take in eqn (5) for dif- 
ferent choices of IMF. Following this method, we find 
that, to match our reference IMF, requires having <j!> — 
1 in equation (5), and multiplying the SFR deduced 
from eqn. (7) by a factor of (1.33/1.02)i "5 = 1.32 to 
account for the fact that Devriendt et al. (1999) use 
a Salpeter IMF with an upper mass limit of 12OM0. 

To calculate the scaling appropriate for the GOfim 
estimator of Cram et al. (1998), we must not only 
account for the different M/Zboi values for their IMF 
and our reference Salpeter (1955) law (given above), 
but also the fact that they quote an SFR which refers 
to stars of mass M > 5Mq, which is done as follows: 
only one ninth of the stellar mass formed with their 
IMF is in stars of M > 5M0 , with the resuh that the 
SFR value from equation (6) must be multiplied by a 
factor 9/(2.52/1.33) — 4.8 to make it appropriate for 
the formation of stars according to our reference IMF. 

The second class of methods comprises the three 
SFR estimators based on monochromatic UV lumi- 
nosities. In this case, the factors of M/Lboi in the 
scaling are replaced by M/Lx, where Lx is the mean 
value of the particular monochromatic luminosity un- 
der consideration, evaluated over the stars formed ac- 
cording to the given IMF. Values for Lx can be readily 
computed under the assumption that main sequence 
stars emit as black bodies in the ultraviolet. This we 
implement, making double and triple power law fits, 
respectively, to the main sequence radius-mass and 
effective temperature-mass relation data tabulated 
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by Binney & Merrifield (1998), i.e.: 



logi 



R 

Rq 

0.02 + 0.72 logi(,(M/MQ) 
0.30 + 0.44 logi(,(Af/MQ) 



M/Mq < 10 
M/Mq > 10 



(A9) 



and 



log 



( 



off 



1" V 103 K 
0.63 + 0.171ogio(M/MQ) 
0.76 + 0.59 logio(M/MQ) 
1.11 +0.30 logio(M/M0) 



(AlO) 

M/Mq < 0.5 

0.5 < M/Mq < 1.6 

M/Mq > 1.6. 



(Note that, if we use these scalings and the rela- 
tion Lboi/i© = (R/R(T,f ( 7;ff/ 5770K)'*. instead of 
the L{M) relation of Eqn. A7, we obtain ratios of 
M/Zboi for pairs of IMF models whic h differ by 
~20%, typically, from those of Eqn. A8.) L\ values 
may then be computed, using the fact that L\ = 
J lx{M)^{M)dM, where lx{M), the monochromatic 
luminosity at wavelength A due to stars of mass M is 
given by Ix = AirRiM)'^ Bx[T^aiM)], where Bx is the 
Planck function, and R and Teg come from Eqns. |A9| 
and IAIOI. 



For the two UV-based estimators of Devriendt et 
al. (1999) we must make conversions at A = 160 and 
280nm, from a Salpeter (1955) IMF with mass range 
[0.1, 120] Af0 to our canonical range of [0.1, 100]Mq. 
The ratio of the M values for the two IMFs is unity 
to better than 1% precision, while the Lx values differ 
by ~ 10%: we obtain scahngs of 1.1^"^ and l.l^ '^^ 
respectively, for the 280nm and 160nm estimators, ac- 
counting for the non-hnear relationship between Lx 
and the SFR from Devriendt et al. (1999). As before, 
we must account for the fact that the Cram et al. 
(1998) U band estimator refers to the formation of 
stars of mass greater than 5Mq only, when we de- 
rive its correction factor, as well as considering the 
difference between the Lx value at 250nm between 
the Cram et al. (1998) IMF and our canonical model. 
M/Lx for the Cram et al. (1998) IMF is only 57 per 
cent of that for the canonical IMF, reflecting their 
difference in slope, so, to correct the SFR estimator 
of eqn (2) we must multiply it by 0.57 x 9 = 5.1. 

The 1.4 GHz estimator of Cram et al. (1998) con- 
stitutes the third class, and it is different in that it 
relates the SFR to the number of stars formed, rather 
than to their luminosity. As presented by Condon 
(1992), the argument for using decimetric radio lumi- 
nosity as a measure of star formation rate hinges on 
the assumption that the radio emission from active 
star-forming galaxies is dominated by non-thermal 
emission from the remnants of Type lb and Type II 
supernovae. In this picture, the radio luminosity is, 
thus, proportional to the (Type lb and Type II) su- 
pernova rate, i^sn, which is equated to the rate of for- 
mation of progenitors of sufficient mass (M > AIq): 



where f^iM > Mq) is the fraction of stars formed 
which have mass M > Mq, and M/M equals the 
number of stars formed per unit time. Since, the deci- 
metric radio luminosity, Lcm, is proportional to i^sN, 
it follows that 



M oc 



M 



fy{M > Mo) 



X Lc 



(A12) 



which gives the correct scaling of M per unit radio 
luminosity between IMFs. If we follow Condon (1992) 
and take Mq = 8Mq, we find that the M/fyiM > 
Mo) value for our canonical IMF is 65 per cent that 
for the IMF assumed by Cram et al. (1998), reflecting 
smaller fraction of high-mass stars in the IMF, 
and from this we find that the factor by which the 
SFR given in eqn (8) has to be multiplied is 0.65 x 9 = 
5.8, where, as before, the factor of nine corrects for 
the fact that eqn (8) only gives the rate of formation 
of stellar mass in stars of mass greater than 5 Mq. 

Similar procedures must be applied to convert 
the SFR estimates presented (for a Salpeter IMF 
over the mass range [l,100]AfQ) by Rigopoulou et al. 
(2000) to our canonical IMF: for their far-infrared es- 
timator, this is an Af/Lboi scaling, while for their Ha 
estimator, we must scale by M /Lx, for A = 6560A. 
Coincidentally, we find that, in both cases, this yields 
a factor of 2.5 by which their SFR values must be mul- 
tiplied for comparison with others computed for our 
canonical IMF. A further correction, however, must 
be applied, since Rigopoulou et al. (2000) assumed a 
cosmology with Q.m =0.3 and Q.a = 0.7, which yields 
a larger luminosity distance out to a given redshift 
than does our canonical Einstein - de Sitter model. 
At the typical redshift {z ~ 0.6), this reduces the SFR 
estimates by ~ 30 per cent, so that the Rigopoulou et 
al. (2000) SFRs are to be multiplied by typical values 
of 1.25 before comparison with our own. 



APPENDIX B: A LOGNORMAL MODEL 
FOR THE SAMPLING BIAS AND 
VARIANCE IN p. 

We want to compute confidence ranges for p*, the 
true, global mean SFR density at z ~ 0.2 and z ~ 0.5 
on the basis of our raw values, given in Section 5.1. 
An accurate determination of these confidence ranges 
requires a model for the galaxy population at < z < 
1 of a sophistication that is beyond the scope of this 
paper. However, we may obtain an estimate for them 
as follows. 

If we denote the SFR density by X, so that the 
desired global mean is X and an estimate is denoted 
by X, then we wish to determine the probability 
distribution p{X\X). We may adopt a Bayesian ap- 
proach, writing 



/^SN = /n(A/ > Mo) X {M/M), 



(All) 



p{X\X) 



p{X,X) _p{X\X)-p(X) 



PiX) 



PiX) 



(Bl) 
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Figure Bl. The probability density functions p(l + 5) as 
a function ofl + (5 = X/X in tlie lognormal model for 
the three cases of =0.5, 1.0 & 2.0, plotted with dotted, 
solid and dashed lines respectively. 



Now, if we assume a uniform prior for p{X) we find 
that the relative probability of X given X is 



p{X\X)^p{X\X)^p{l + S), 



(B2) 



where 5 is the fractional den sity fluctuation. Since our 
estimates from Section 5.2 for the variances in our 



two sampling volumes are of order unity, we cannot 
assume that the distribution of 5 is Gaussian. Indeed, 
gravity will have skewed the distribution so that it is 
significantly more likely that a randomly-located sur- 
vey volume will be underdense than overdense, so we 
are more likely to have under-estimated p* than to 
have over-estimated it. For present purposes, we as- 
sume a lognormal model for the density field, because 
Bernardeau & Kofman (1995) argue that it gives a 
good fit to the probability density function of the 
density fluctuations in TV-body simulations of a Cold 
Dark Matter universe well into the non-linear regime 
(i.e. at least until the variance in fractional overden- 
sity reaches ~ 1.5^). 

In the lognormal model, we have 



p(l + 5)d(l + 5) 



1 



X exp 



(In(l + 5) + Gl/2f 



2ol 



A{l + 5) 

+ ' 



(B3) 



where ctq = ln(l -I- cr^) is the variance of the ficti- 
tious Gaussian field from which the true density field, 
with variance cr^, is supposed to be derived, by tak- 
ing the natural logarithm. It can readily be shown 
that p{l -\- 5) takes its maximum value at 1 -I- 5 = 
exp(— 3(To/2), where p{l + 5) = exp(ao ) / (27rCTQ ) . 



In Fig. gl| we plot p{l + 5) against 1 + 5 = X/X 
for the three cases of cr^ = 0.5 (dotted line), 1.0 (solid 
line) and 2.0 (dashed line), assuming the lognormal 
model of eqn. B3. The relatively small fraction of 
area under the curves at 1 + 5 > 1, indicates the 



low probability that the true global mean SFR den- 
sity will be lower than that estimated, reflecting the 
fact that for > 0.5 much more of the universe is 
underdense than overdense. However, the mean value 
oi 1 + 5 must always equal unity, by definition, so our 
estimated SFR density is not a biased estimate, but 
the evolution of the cosmological density field through 
gravitational instability does act to increase the width 
of the confidence range for p* at a given percentile. 

For the lognormal model, this is simple is to com- 
pute. To obtain the width of, say, the 68 per cent 
confidence region for 1 -I- 5, we want to take the dif- 
ference Z2 — zi where z = 1 + 5 — zi and z = Z2 are 
the solutions to the equation 



p{l + 5)d(l + 5) =p 



(B4) 



for p=0.16 and 0.84 respectively. It is easy to show 
that 



p{l+5)d{l+5) = - 



1 -l-erf 



hi{z) + al/2 



,(B5) 



where erf(2;) is the error function of x, and, hence, 
to obtain 68 per cent confidence 
ranges of 0.65 <X/X < 2.3, 0.62 <X/X < 3.2, and 
0.61 <X/X < 5.0, for ^^=0.5, 1.0 and 2.0, respec- 
tively. It is interesting to note that, in this lognor- 
mal model, at least, the effect of continued evolu- 
tion of the density field via gravitational instability 
is really only seen in the shift of the upper limit to 
the confidence interval for X/X, and that the lower 
limit barely changes: so, as time passes and in- 
creases, the amount by which one's survey may have 
under-estimated the true, global mean SFR density 
increases markedly, but the amount by which it may 
have over-estimated it stays pretty much the same. 
This is equivalent to the noting that the difference 
between the three curves in Fig. 



Bl 



is much more 
pronounced at low values oi 1 + 5, while, at high val- 
ues, they are very close. 
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